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Dedicated to Bill Thurston on his 60th birthday. 

Abstract. In this paper we develop a new theory of infinitesimal harmonic 
deformations for compact hyperbolic 3-manifolds with "tubular boundary". 
In particular, this applies to complements of tubes of radius at least Ro = 
arctanh(l/v / 3) ~ 0.65848 around the singular set of hyperbolic cone manifolds, 
removing the previous restrictions on cone angles. 

We then apply this to obtain a new quantitative version of Thurston's 
hyperbolic Dehn surgery theorem, showing that all generalized Dchn surgery 
coefficients outside a disc of "uniform" size yield hyperbolic structures. Here 
the size of a surgery coefficient is measured using the Euclidean metric on 
a horospherical cross section to a cusp in the complete hyperbolic metric, 
rescaled to have area 1. We also obtain good estimates on the change in 
geometry (e.g. volumes and core geodesic lengths) during hyperbolic Dehn 
filling. 

This new harmonic deformation theory has also been used by Bromberg 
and his coworkers in their proofs of the Bers Density Conjecture for Kleinian 
groups. 



1. Introduction 

Let X be a compact, orientable 3-manifold with a finite number of torus bound- 
ary components and suppose its interior admits a complete, finite volume hyperbolic 
metric. For each torus, there are an infinite number of topologically distinct ways to 
attach a solid torus, corresponding to the homotopy class of the non-trivial simple 
closed curve on the boundary torus that bounds a disk in the solid torus. The set 
of non-trivial simple closed curves on a torus is parametrized by pairs of relatively 
prime integers, once a basis for the fundamental group of the torus is chosen. If 
each torus is filled, the resulting manifold is closed. A fundamental theorem of 
Thurston ([18]) states that, for all but a finite number of filling curves on each 
boundary component, the resulting closed 3-manifold has a hyperbolic structure. 
Those that don't result in a hyperbolic structure are called exceptional curves. 
In we showed that there is a universal bound to the number of curves that 
must be excluded from each boundary component; in particular, when there is one 
boundary component, there are at most 60 exceptional curves and when there are 
multiple boundary components, at most 114 curves from each component need to 
be excluded. 

The proof of this result uses the harmonic deformation theory developed in [10 
to deform the finite volume complete hyperbolic structure on the interior of X. 
The deformation consists of a family of singular hyperbolic metrics on the filled 
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manifolds where the singularities lie along geodesies isotopic to the cores of the 
attached solid tori. The metrics on discs perpendicular to the geodesies have a 
single cone point at the geodesic with a cone angle a which is constant along each 
component. We call such structures hyperbolic cone manifolds. At the complete 
structure the cone angles are considered to be equal to and if the deformation 
reaches cone angles 2n on each component, this is a smooth hyperbolic structure 
on the filled manifold. 

For the analytic techniques in [TU] to work, it is necessary to restrict all the cone 
angles to be at most 2ir. While this is adequate for the Dchn filling problem, there 
are other situations where it is important to be able to deal with larger cone angles 
or with a more general type of singularity. In particular, in his proof of the Bers 
Density Conjecture, Bromberg ([4], [5]) (and then Brock-Bromberg ([6]) for more 
general versions of the Density Conjecture) needs to deform cone angles equal to 47T 
back to 2-7T. In order to have a similar deformation theory in this type of situation, 
a new analytic technique is necessary. One of the primary goals of this paper is 
to provide such a technique. To this end, the previous analysis, which developed 
an L 2 Hodge theory for the incomplete smooth metric on the complement of the 
singular locus, is replaced by a Hodge theory on the compact hyperbolic manifold 
with boundary obtained by removing an open tubular neighborhood of the singular 
locus. 

We will now give a brief outline of this theory and some of its applications, as 
they are presented in the body of this paper. 

If one removes an equidistant tubular neighborhood of each component of the 
singular locus of a 3-dimensional cone manifold, one obtains a smooth hyperbolic 
manifold with torus boundary components. The boundary tori have intrinsic flat 
metrics. Furthermore, the principal normal curvatures are constant on each com- 
ponent, equal to n, ^ (where we assume that k > 1). If n ^ 1 the lines of curvature 
are geodesies in the flat metric corresponding to the meridional and longitudinal di- 
rections. The normal curvatures and the tube radius, R, are related by coth R = k 
so they determine each other. We call tori with these curvature properties tubular. 

We say that an orientable hyperbolic 3-manifold M has tubular boundary if its 
boundary components are flat tori with constant principal normal curvatures as 
above. The interior of M is known to also admit a complete finite volume hyper- 
bolic metric ([TTJ Lemma 3.8]). Associated to M is a filled hyperbolic manifold 
M. It is canonically obtained from M by extending the hyperbolic structures on 
the boundary tori normally as far as possible. Each added component is topologi- 
cally a torus crossed with R and is foliated by tubular tori whose radii go to zero. 
Then M is homeomorphic to the interior of AI and has a (typically incomplete) 
hyperbolic structure. The original manifold with tubular boundary M is a subset 
of M obtained by truncating the ends of M along the appropriate tubular tori. 
Note that M contains many other manifolds with tubular boundary obtained by 
choosing other truncations. 

The geometry of the ends of the filled manifold M is completely determined by 
the geometry of the tubular boundary tori of M. In particular, if, for a particular 
boundary component, the lines of curvature for the larger curvature k > 1 are 
parallel closed curves, the completion of that end of M will have the structure of a 
hyperbolic cone manifold with those curves as meridians around the singular locus. 
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The resulting cone angles can be read off from the geometry of the boundary and 
can be arbitrarily large. In general, these lines of curvature merely determine a 
foliation on the boundary torus where the leaves are geodesic in the flat metric on 
the torus. There is still a canonical way to extend the structure of the boundary 
torus in this case but the singular set for the completion will be a single point with 
a complicated neighborhood. The resulting structure obtained by completing M, 
including both cone manifolds and this more general type of singularity, is called a 
hyperbolic structure with Dehn surgery singularities. (See |18j for details.) 

Fix a component T of the boundary of a hyperbolic 3-manifold M with tubu- 
lar boundary and consider the holonomy group of the fundamental group of that 
boundary torus. Assume that n > 1; then each element in the holonomy group will 
have an invariant axis. Since the group is abelian, all of the elements fix a common 
axis in H 3 . Choose a direction along the axis. Then, associated to each element is 
a complex length whose real part is the signed translation distance along the axis 
and whose imaginary part is the total rotation around the axis. The amount of 
rotation is a well-defined real number whose sign is determined by the right hand 
rule. The map sending an element of ni(T) — Hi(T;T,) to its complex length is 
linear and can be extended canonically to a linear map £ : Hi(T;M) — » C. The 
resulting value C(c) for any element c £ Hi(T; R) will be called the complex length 
of c. 

When the complex length of a simple closed curve 7 on T equals 27ri, this implies 
that the end of the filled manifold M corresponding to T completes to a smooth 
structure on the manifold M (7) obtained from M by Dehn filling with filling curve 
7. More generally, if the complex length of 7 equals m, the end completes to a 
hyperbolic cone structure on M(pf) with cone angle a. Now assume that the map 
£ : Hi(T;M) — > C is a (real) isomorphism. (This holds whenever the complex 
lengths of two generators of f/i(T;Z) are linearly independent over the reals, a 
condition that holds whenever M has tubular boundary with k > 1 on T.) Then 
there will be a unique element c £ H\{T\ R) so that £(c) = 2wi; we say that c is the 
Dehn surgery coefficient of the boundary component T. When c £ Hi(T;Z) is a 
primitive element, it corresponds to a simple closed curve and the completion of the 
end is smooth. When it is in Hi(T;Q), the completed end has a cone singularity 
along a core geodesic. 

Suppose, for simplicity, that M has a single boundary component T and denote 
by X the underlying smooth manifold with boundary. Then the subset of Hi (T; R) 
consisting of Dehn surgery coefficients of hyperbolic structures with tubular bound- 
ary on X is called the hyperbolic Dehn surgery space for X and will be denoted 
by HT)S(X). (Note: In [IS] hyperbolic Dehn surgery space is equivalently defined 
in terms of the filled in structures.) Thurston's theorem about the finiteness of 
exceptional curves is actually a corollary of his theorem that HT)S(X) contains a 
neighborhood of infinity (infinity here corresponds to the complete finite volume 
structure on the interior of X). Hence, it contains all but a finite number of points 
of the integral lattice Hi(T;Z) c i?i(T;R). Similar statements are proved when 
there are multiple boundary components. 

Thurston's proof is not effective; it gives no information about the size or shape 
of hyperbolic Dehn surgery space which is why there is no information from his 
proof about the number of exceptional curves in his finiteness theorem. Note that 
the vast majority of the points with integral entries are "near" infinity and thus 
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that the statement that HT>S(X) contains a neighborhood of infinity is really a 
much stronger statement than the finiteness of exceptional fillings. 

One of the main goals of this paper is to provide an effective proof of Thurston's 
result, one that will guarantee that TCDS(X) contains a neighborhood of infinity 
of a "uniform" size and shape. To simplify the description of this uniform region, 
it is useful to put a metric on Hi(T;M.). One way to do this is to consider the 
complete structure on the interior of X and take a horospherical torus T embedded 
in its end. This torus inherits a flat metric which is well-defined up to scale. The 
homology group H±(T;M.) can be canonically identified with the universal cover of 
T. So the flat metric on T, normalized to have area 1, induces a flat metric on 
Hi(T;M.) = R 2 . Note that under this identification, the distance from the origin 
to point c € Hi(T; Z) is just the geodesic length of the corresponding closed curve 
7 on T measured with respect to the flat metric on T, normalized to have area 1. 
This is called the normalized length of 7 on T. As with complex length this notion 
of length can be extended naturally to define a map L : Hi (T; R) — > R; the value 
I/(c),c G i?i(T;R) is called the normalized length of c. 

The theorem below says that, using this metric on the plane, TCDS(X) always 
contains the complement of a disk of uniform radius around the origin, independent 
of X. 

Theorem 1.1. Consider a complete, finite volume hyperbolic structure on the in- 
terior of a compact, orientable 3-manifold X with one torus boundary component. 
Let T be a horospherical torus which is embedded as a cross-section to the cusp of 
the complete structure. Consider TU)S(X) as a subset of Hi(T;M) = R 2 where 
the latter is endowed with the Euclidean metric induced from the universal cover 
of T with its flat metric scaled to have unit area. Then TCDS{X) contains the 
complement of a disk of radius 7.5832, centered at the origin. Equivalently, any 
c 6 Hi(T;M) whose normalized length L{c) is bigger than 7.5832 is in TCDS(X). 

In [11] we showed that any simple closed curve 7 on T, viewed as an element of 
-ffi(T;Z), whose normalized length is at least 7.515 is in 7iT>S(M). Thus, except 
for the slight change in constant (which is due to the tube radius condition required 
for the Hodge theorem, as discussed below), Theorem 1 1.1 1 is a direct generalization 
of that result. The normalized length condition translates easily into an upper 
bound on the number of exceptional fillings. 

The proof of Theorem 11.11 like the proof of the uniform bound on exceptional 
fillings, involves two main steps. First, it is necessary to show that one can deform 
a given structure towards the desired structure. For example, in the cone manifold 
case, one needs to show that the cone angles can always be increased a small 
amount. This step depends on proving a local parametrization theorem, showing 
that, locally, the deformations are parametrized by their Dehn surgery coefficients. 
In order to find such a local parametrization, one proves a local rigidity theorem 
which says that it is impossible to deform the hyperbolic structure while keeping 
the Dehn surgery coefficients fixed. The local parametrization then follows by an 
application of the implicit function theorem. 

The second step is to show that, under certain initial conditions, it is always 
possible to deform the complete structure on the interior of X, through hyperbolic 
structures with tubular boundary, to one with the desired Dehn surgery coefficient 
before there is any degeneration of the hyperbolic structure. This step requires 
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one to control the change in geometry under the deformation and depends on the 
analysis in the proof of local rigidity. 

The proofs of the local rigidity and local parametrization theorems require new 
analytic techniques and occupy the next three sections. Once these are established, 
the arguments to establish uniform bounds closely follow those in However, 
the use of manifolds with tubular boundary, as opposed to cone manifolds, leads to 
subtly different estimates when there are multiple boundary components. 

For a compact, orientable, 3-manifold X with multiple torus boundary compo- 
nents Ti, • ■ • , Tfe, the Dehn surgery space is a subset of ®iHi(Ti\M.) and the Dehn 
surgery coefficient c = (ci,C2,--- , c&) is determined by the Dehn surgery coeffi- 
cients Cj £ Hi(Ti;R) for each torus. For ci £ Hx(Ti;M) the normalized length L(ci) 
is computed on a horospherical torus corresponding to Ti in the complete structure 
on the interior of X as described above. In this case we prove the following uniform 
statement: 

Theorem 1.2. Consider a complete, finite volume hyperbolic structure on the inte- 
rior of a compact, orientable 3-manifold X with k > 1 torus boundary components. 
Let Ti,-- - ,Tfe be horospherical tori which are embedded as cross-sections to the 
cusps of the complete structure. Consider TCDS{X) as a subset of ©jif^T^R). 
Then there exists a universal constant C = 7.5832 such that c = (ci, C2, • • • , c&) is 
in TCDS{X) provided the normalized lengths Li = L(ci) satisfy 



When k = 1 this is precisely the same statement as that of Theorem 11.11 In 
the multiple cusp case, it again gives a uniform upper bound on the number of 
exceptional simple closed curves that need to be excluded from each boundary 
component so that the remaining Dehn filled manifolds are necessarily hyperbolic. 

However, it should be noted that the bound depends on the number of boundary 
components. This is in contrast to Theorem 5.12 in [11] which provides a uniform 
bound independent of the number of cusps. The reason for this difference is that in 
the previous paper we allow the possibility of increasing the cone angles at varying 
rates. Once a cone angle of 2ir is attained on one component of the singular locus, 
it is no longer changed, while the other angles are increased. We no longer keep 
track of the geometry in a neighborhood of the smooth core geodesies. (Indeed, 
a geodesic could become non-simple and change isotopy class.) This was not ade- 
quately explained in [llj ; for a discussion of this and other subtler issues that arise 
in the multiple cusp case, the reader may consult [IB] . 

In the current paper, we require a lower bound on the tube radius of all the 
tubular boundary components throughout the deformation. This is because, even at 
the final Dehn surgery coefficient, the filled manifold may still have singularities, and 
thus we no longer have the luxury of ignoring the tube radius around a component 
once its desired surgery coefficient is attained. We always move radially in HT>S(X) 
from the complete structure to the desired Dehn surgery coefficients (which, when 
all the coefficients correspond to simple closed curves, amounts to increasing the 
cone angles at equal rates) . This provides weaker estimates in the case of multiple 
cusps. It does, however, have the advantage in the case of smooth Dehn filling 
(or cone manifolds) that the isotopy class of the union of the core geodesies will 
necessarily remain unchanged. 
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We also obtain good control on the change in geometry during generalized hy- 
perbolic Dehn filling. Theorem 15.121 gives explicit upper and lower bounds for the 
volume and core geodesic length, with the asymptotic behavior given by Neumann- 
Zagier in [15] . These bounds are illustrated in Figures [2] and [3] at the end of the 
paper. For example, we obtain the following numerical bound. 

Theorem 1.3. Let X be a compact, orientable 3-manifold as in Theorem \ 1 . H and 
let Vac denote the volume of the complete hyperbolic structure on the interior of X . 
Let c = (ci, . . . ,Cfe) € Hi(dX;M) be a surgery coefficient with normalized lengths 
Li = L(c.i) satisfying 

V — < 4r where C = 7.5832, 
L 2 C 2 

and let M(c) be the filled hyperbolic manifold with Dehn surgery coefficient c. Then 
the decrease in volume AV = — vol(M(c)) during hyperbolic Dehn filling is at 
most 0.198. 

We now briefly explain how the use of manifolds with tubular boundary allows 
us to avoid the analytic issues that led to the cone angle restriction in [IT] . 

The original local rigidity theory in [TO] applies only to hyperbolic 3-manifolds 
with conical singularities along a geodesic link where the cone angles are restricted 
to be at most 2ir. The argument involves finding, for any infinitesimal deformation 
of the hyperbolic cone structure, a harmonic representative and then utilizing a 
Weitzenbock formula for such harmonic infinitesimal deformations. The analysis 
using this formula involves an integration by parts on the complement of a tubular 
neighborhood of the singular locus, resulting in a term from the boundary of the 
tubular neighborhood. Any deformation for which this boundary term goes to zero 
as the radius of the tube goes to zero is seen to be trivial. The main step in [TU] is 
to show that, for any infinitesimal deformation where all the cone angles are held 
constant, the boundary term does, indeed, go to zero as long as the cone angles are 
at most 2ir. The fact that the analysis involves arbitrarily small neighborhoods of 
the singular locus means that it depends on the asymptotic behavior of harmonic 
deformations near the singular locus. This behavior is strongly governed by the 
value of the cone angle along the singular locus. 

To avoid a dependence on the local behavior near the singular locus, it is nec- 
essary to work on the complement of a tubular neighborhood whose tube radius is 
bounded below. Again, one must find a harmonic representative for the infinites- 
imal deformation of the hyperbolic structure with tubular boundary. Since these 
infinitesimal deformations can be viewed as cohomology classes, this can be viewed 
as Hodge theory on a manifold with boundary. In this case, when one uses the 
Weitzenbock formula and integration by parts, one wants to end up with a bound- 
ary term with an appropriate sign. When such a boundary term is obtained, the 
conclusion is again that the deformation is trivial, implying a local rigidity theorem 
as before. 

Thus, a Hodge theory must be developed with the boundary term from the 
Weitzenbock formula in mind. To this end, we find a formula for this term in 
Section[3] The formula derived there is quite general and is valid for any hyperbolic 
3-manifold with boundary, not just those with tubular boundary. As a result, it 
should have applications in other contexts and may be of independent interest. 
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The form of this Weitzenbock boundary term motivates the specific boundary 
conditions we require for our Hodge representative when the manifold has tubular 
boundary. The proof that the corresponding boundary value problem can always 
be uniquely solved is contained in Section [4j This result requires a universal lower 
bound on the tube radius of the tubular boundary components. By definition 
the tube radius R of a tubular boundary component is determined by the formula 
coth R — k, where k > 1 is the larger of the two principal curvatures on that 
component. When k — 1 (which corresponds to a horosphcrical torus), the tube 
radius is said to be infinite. 

Once the required Hodge theorem is proved, similar arguments to those in [10] 
imply the following local rigidity and local parametrization result. The previous 
cone angle restriction has been removed and is replaced by a mild restriction on 
the tube radius. For simplicity, we also assume that all tube radii are finite. 

Theorem 1.4. Let M be a compact, orientable hyperbolic 3-manifold with tubular 
boundary and suppose that the tube radius of each boundary component is finite 
and at least Rq — arctanh(l/v / 3) ~ 0.65848. Then there are no deformations of 
the hyperbolic structure fixing the Dehn surgery coefficient of M . Furthermore, 
the nearby hyperbolic structures with tubular boundary are parametrized by their 
Dehn surgery coefficients. In particular, a finite volume hyperbolic cone-manifold 
with singularities along a link and tube radii at least Ro has no deformations of 
the hyperbolic structure keeping the cone angles fixed, and the nearby hyperbolic 
cone-manifold structures are parametrized by their cone angles. 

Once this analytic theory for hyperbolic manifolds with tubular boundary is 
developed and the above local rigidity theorem is proved, the arguments in go 
through with minor changes. Indeed, much of that paper was written in the context 
of manifolds with tubular boundary, once the necessary analytic and geometric 
control was derived. These arguments are recalled in Section where they are 
then applied to prove Theorems 11.11 and 11.21 and other results. 

2. Preliminary Material 

In this section we recall the basic setup for the harmonic deformation theory of 
hyperbolic structures on 3-manifolds. The reader is referred to the papers [10] and 
[iT] for details, and to [T^] for a survey of the theory and its applications. 

An infinitesimal deformation of a hyperbolic structure on a hyperbolic 3-manifold 
M is given by a cohomology class in H 1 (M\ E) where is E is the bundle of (germs 
of) infinitesimal isometries of H 3 . By viewing this cohomology group in terms of 
de Rham cohomology, such a cohomology class can be represented by a 1-form 
with values in E. The 1-form will be closed with respect to the E'-valued exterior 
derivative which we denote by dE- A representative for a cohomology class can 
be altered by a coboundary without changing its cohomology class. An _E-valued 
1-form is a coboundary precisely when it can be expressed as dES, where s is an 
E- valued 0-form, i.e. a global section of E. 

A standard method for choosing a particularly nice representative in a coho- 
mology class is to find a harmonic representative: one that is co-closed as well as 
closed. On a closed manifold such a harmonic representative is unique. When the 
manifold is non-compact or has boundary, it is necessary to choose asymptotic or 
boundary conditions to guarantee existence and uniqueness. 
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When M is a hyperbolic 3-manifold with tubular boundary, one can begin with 
a representative Co £ H l {M]E) that has a special form in a neighborhood of the 
boundary. In [10) , when M is a hyperbolic cone manifold, specific closed E- valued 
1-forms, which we call standard forms, are defined in a neighborhood of the singular 
locus. The same forms are defined in the neighborhood of the boundary components 
of a general hyperbolic manifold with tubular boundary. They have the property 
that any possible infinitesimal change in the holonomy representation of the fun- 
damental group of a boundary torus can be induced by one of these forms. As 
a result, by standard cohomology theory, for any infinitesimal deformation of the 
hyperbolic cone manifold structure, it is possible to find a closed -E-valued 1-form 
Co on M which equals one of these standard forms in a neighborhood of each torus 
boundary. 

The standard forms are harmonic so the _E-valued 1-form Co will be harmonic in 
a neighborhood of the boundary but not generally harmonic on all of M. Since it 
represents a cohomology class in H 1 (M; E), it will be closed as an E- valued 1-form, 
but it won't generally be co-closed. If we denote by 6e the adjoint of the exterior 
derivative, d,E, on E- valued forms, then this means that d,E& = 0, but SeCo 7^ 
in general. Finding a harmonic (i.e. dE closed and Se co-closed) representative 
cohomologous to Co is equivalent to finding a global section s such that 

(2) S E d E s = S E Co. 

Then, lo = Co + <igs satisfies Se^> = 0, d^a; = 0; so it is a closed and co-closed 
representative in the same cohomology class as Co. 

The fibers of the bundle E are all isomorphic to the Lie algebra g of the Lie group 
of isometries of hyperbolic space. One special feature of the 3-dimensional case is 
the complex structure on the Lie algebra q = sfaC The infinitesimal rotations 
fixing a point p £ H 3 can be identified with su(2) = so(3), and the infinitesimal 
pure translations at p correspond to i su(2) = T p H 3 . Geometrically, if t £ T p M 3 
represents an infinitesimal translation, then i t represents an infinitesimal rotation 
with axis in the direction of t. Thus, on a hyperbolic 3-manifold M we can iden- 
tify the bundle E with the complexified tangent bundle TM (g> C. At each point, 
the fiber decomposes into a real and imaginary part, representing an infinitesimal 
translation and an infinitesimal rotation, respectively, and we can speak of the real 
and imaginary parts of an _B-valued form. 

In [10] it was shown that in order to solve equation for i?-valued sections, 
it suffices to solve the real part of the equation. The real part of a section s of E 
is just a (real) section of the tangent bundle of M; i.e., it is a vector field, which 
we denote by u. The real part of 8e<1es equals (V*V + 2) u, where V denotes the 
(Riemannian) covariant derivative and V* is its adjoint. The composition V*V is 
sometimes called the "rough Laplacian" or the "connection Laplacian" . 

To solve the real part of the equation ((2J), we find that the computations are 
somewhat easier if we replace vector fields by their dual real-valued 1-forms. We 
take the real part of —SeCo, considered as a vector field, and denote its dual 1-form 
by C- The operator (V*V+2) on vector fields becomes A + 4 on dual 1-forms, where 
A = dS + 6d is the usual Laplacian on real- valued 1-forms and we are denoting the 
exterior derivative and its adjoint on M by d and S, respectively. We must then 
solve the equation 



(3) 



(A + 4)r = C, 
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for a globally defined real- valued 1-form ronM, which will be dual to the vector 
field u. 

For a manifold with boundary it is necessary to prescribe boundary conditions 
on r for this problem to have a unique solution. The boundary conditions we choose 
are non-standard and very specific to our hyperbolic deformation theory context. In 
particular, the local rigidity results that we seek depend on a Weitzenbock formula 
for harmonic -E-valued 1-forms. This formula contains a boundary term whose 
sign is crucial to the argument. The harmonic form lo G H (M; E) is obtained 
by solving ([3]), which, in turn, gives us a solution to Since lo = lo + cIeS and 
lo is standard in a neighborhood of the boundary, our boundary term will have a 
contribution from the standard form, which is quite explicit, and from the correction 
term cIes. The behavior of the latter depends on our choice of boundary condition 
when solving (|3|). A major consideration when choosing a boundary condition is 
that the contribution to the Weitzenbock boundary term from the correction term 
cIes be non-positive. 

In Section[3]we compute a formula for this contribution for a general hyperbolic 3- 
manifold with boundary. We then specialize to our current situation of a hyperbolic 
manifold with tubular boundary and choose boundary conditions specific to this 
case. In Section 2] we prove that the problem of solving © with these boundary 
conditions always has a unique solution. 

We now recall the Weitzenbock formula for harmonic E- valued 1-forms, referring 
to [10] and [UJ for details and proofs. 

We can decompose any lo € H 1 (M,E) into its real and imaginary parts lo = 
rj + if), where rj and rj are vector field valued 1-forms on M which we can view 
as elements of Horn (TM, TM) at each point of M. The real symmetric part of 
lo, viewed as a symmetric 2-tensor, describes the infinitesimal change in the metric 
induced by the infinitesimal deformation corresponding to lo. 

One can always choose a representative for a cohomology class where rj is sym- 
metric, when viewed as a section of the bundle Horn (TM, TM ). To do this, one 
notes that, since lo is d^-closed, it is the image of a locally defined section. Then rj 
is symmetric if the local section has the property that its imaginary part equals | 
of the curl of the real part, where both the real and imaginary parts are viewed as 
locally defined vector fields. It is shown in [10] that such a choice of local section 
is always possible; in that paper, such a local section was called a canonical lift of 
the real part. (Note, however, that the definition of curl in that paper differs from 
the standard one, which is the one used here, by a sign and a factor of 2). 

If lo is harmonic then rj satisfies the equation 



where D denotes the exterior covariant derivative on vector valued 1-forms and D* 
is its adjoint. If rj is also traceless then it satisfies D*rj = 0; hence, it satisfies the 
simpler equation 



In this case we have that rj = iDrj, where i is the Hodge star operator on forms 
in M and takes the vector valued 2-form Dr\ to a vector valued 1-form. It is also 
true in this case that *Drj is symmetric and traceless. Thus, we can write 



(D*D + DD*)i) + r/ = 



(4) 



D*Drj + n = 0. 



(5) 



LO = 1] + f) 



rj + i iDrj, 
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where both n and *Dw are symmetric and traceless. 

Let M be a hyperbolic 3-manifold with tubular boundary, whose boundary com- 
ponents are tori of tube radii R\ , . . . , Rk ■ The boundary components will always be 
oriented by the inward normal for M. For any TM-valued 1-forms a, f3 we define 



where R denotes the vector (i?i, . . . , Rk). 

In this integral, *Da A j3 denotes the real valued 2-form obtained using the 
wedge product of the form parts, and the geometrically defined inner product on 
the vector- valued parts of the TM-valued 1-forms *Da and (3. 

Returning to equation |4]), we take the L 2 inner product on M of this equation 
with n and integrate by parts. We then obtain the following Weitzenbock formula 
with boundary for any harmonic infinitesimal deformation lu of the form ([5|): 



In particular, for a non-trivial infinitesimal deformation, the boundary term 
bn(r],r]) must be positive. The proof in [lOj that there are no infinitesimal de- 
formations of hyperbolic cone manifolds (with cone angles at most 2n) fixing the 
cone angles amounts to showing that, for a deformation fixing the cone angles, a 
harmonic representative ui can be found so that this boundary term goes to as 
the tube radius goes to zero. The results in the current paper depend on showing 
that when the Dehn surgery coefficients are all preserved (and the tube radii are all 
bigger than a universal constant), a harmonic representative can be found so that 
this boundary term is non-positive. 

Recall that the harmonic form lu will be found by starting with a representative 
lu which equals a standard form luq in a neighborhood of the boundary, hence is 
harmonic in that neighborhood, but not globally harmonic. Solving equation ([2|) 
provides a correction term cLes that is added to lu to make it globally harmonic. We 
denote the correction term in a neighborhood of the boundary by lu c and decompose 
the harmonic representative lu as lu = luq+lu c in that neighborhood. Note that, since 
luq is harmonic, the correction term lu c will also be harmonic in that neighborhood. 

The standard forms are all of the form (J5j) and their real parts satisfy equation 
(0]). Thus, being able to write lu in this form is equivalent to being able to solve @ 
in such a way that lu c can be written in this form. Then, in a neighborhood of the 
boundary, we write the real part of lu as the sum of the real part of luq and that of 
lu c , 7] — rjQ + n c . Both 7Jq and r\ c satisfy equation |4]) in that neighborhood. 

Using this decomposition of r\ on the boundary, we can try to compute the 
boundary term 6^(77,77). In [TT] we saw that the cross-terms vanish so that the 
boundary term is simply the sum of two boundary terms: 



(6) 




(7) 



\\DvW 2 + INI 2 



(8) 



bR{v,V) = bR{q ,ri ) + b R {rj c ,rj c ). 



Since the standard forms are quite explicit, it is fairly easy to find conditions 
under which the term bu(j]o, 770) is non-positive and to estimate its value in general. 
Thus, we finally come to the boundary value problem we wish to solve: 
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Boundary Value Problem: Find boundary conditions on the real-valued 1-form 
t so that there is always a unique solution to equation (0) when £ is smooth on all of 
M , including the boundary. Furthermore, these boundary conditions must ensure 
that tu c satisfies {5p, hence that rj c satisfies equation Q). Finally, the boundary 
term, bn(ji c ,r] c ), in the Weitzenbock formula must always be non-positive. 

3. The Weitzenbock correction term 

In this section we derive a general formula for a boundary integral (see (| 1 1|) 
below) which we refer to as the Weitzenbock correction term, that arises in the 
Weitzenbock formula ([7]) for harmonic infinitesimal deformations of a compact hy- 
perbolic 3-manifold with boundary. We then specialize to the special case of interest 
in this paper, when the boundary is tubular. We further compute the boundary 
term in this case under the hypothesis of specific boundary conditions. In the next 
section we show that such boundary conditions can always be realized. 

Let M be an oriented compact hyperbolic 3-manifold with boundary dM and 
let E be the bundle of (germs of) infinitesimal isometries on M . We denote by 
dE the coboundary operator on smooth E- valued i-forms on M; the latter are 
denoted by fi*(M; E). This operator satisfies the equation d% = and H 1 (M; E), 
the first cohomology of M with coefficients in E, is defined to be the ds-closed E- 
valued 1-forms modulo those of the form dES where s is an E- valued 0-form; i.e., a 
global section of E. This cohomology group represents the (scheme of) infinitesimal 
hyperbolic deformations of M. 

As discussed in the previous section, a boundary integral, bnfn, rf), occurs in the 
Weitzenbock formula that holds for a class of harmonic (d^-closed and co-closed) 
E- valued 1-forms. In the case when M has tubular boundary, these harmonic 
forms are constructed by adding a coboundary of the form dES to a representative 
in H 1 (M; E) which is in a standard form near the boundary. In particular, we are 
interested in the contribution to the boundary integral coming from this cobound- 
ary. Because of the decomposition ([8]) of the boundary integral, this contribution 
can be computed as a boundary integral involving only the E- valued 1-form d^s. 
In this section we compute this boundary integral, on a general compact hyperbolic 
3-manifold with boundary, for any E- valued 1-form that is of the form dES. 

In the previous section we observed that the bundle E can be identified with 
the complexified tangent bundle TM £g> C. Then, a global section s of E can be 
written as s = u + ui where u,u are global vector fields on M. Similarly we can 
decompose rf^s € CI (M ; E) into its real and imaginary parts d^s = rj c + f/ c i. Both 
f] c and ij c are vector field valued 1-forms; i.e., elements of fi 1 (M;TM). They can 
equivalently be viewed as elements of Horn (TM,TM). In [10] we computed that 

(9) dES = n c + rj c i = {Du — rot^) + (Du + vot u )i 

where Du £ il 1 (M; TM) is the covariant derivative of the vector field u and rot u £ 
Horn (TM, TM) at any point p 6 M is the infinitesimal rotation determined by the 
tangent vector u(p) € T p M. Thus for each tangent vector X £ TM, 

(10) Du(X) = V x u and rot u (X) = u x X, 

where x denotes the cross product defined by the Riemannian metric and orienta- 
tion on M. 
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We now define the boundary integral of interest to us. For any element cIes — 
f]c + fjd £ f2 1 (M; E) we define the Weitzenbock correction term by 

(U) b = [ Vc Af] c , 

JdM 

where the boundary is oriented with respect to the inward normal. Of course, this 
integral can be defined for any element of Q}(M\ E), decomposed into its real and 
imaginary part. However, we will only be interested in this section in computing 
it for those elements which are coboundaries; hence, the name "correction term" . 
Much of our computation is valid for any such element, but we will then specialize 
to the case where both r\ c and fj c are symmetric which is the case that will arise 
during the process of finding a harmonic representative discussed in the previous 
section. In that situation we will also have the relation fj c = iDrj c from which 
it follows immediately that the Weitzenbock correction term equals — bn(r) c ,r] c ), 
where 6^ is defined by ([6]). As discussed in the previous section, because of the 
Weitzenbock formula (JJJ and its decomposition (|8|) we will be interested in finding 
boundary conditions on s that will guarantee that bn(r] c ,r] c ) is non-positive. Thus, 
we will be interested in conditions that will imply that the Weitzenbock correction 
term is non-negative. 

In order to compute this boundary integral, it is useful to decompose sections 
and forms into their tangential and normal parts near the boundary. Specifically, 
the surfaces equidistant from dM give a foliation in a neighborhood of dM in M, 
and there is a unit vector field n consisting of normal vectors to these equidistant 
surfaces pointing inwards from dM. The vector field u on M can be decomposed 
near dM as 

u = v + hn 

where v is the component tangent to the equidistant surfaces and h = u ■ n is the 
component in direction of the normal n. (We use • to denote the Riemannian inner 
product on M .) Similarly, we write 

u = v + hn. 

As in ([9]) above we write oIeS = r/ c + fj c i where s = u + ui is a global section, 
decomposing both s and oIeS into their real and imaginary parts. We are only 
interested in the values of r\ c and fj c restricted to dM. Viewed as a TM-valued 
1-form there, for each X € T p (dM) we decompose n c (X) into a tangential part 
G(X) and normal part F(X). We can then write 

n c (X) = G(X) + F(X)n, where G(X) S T p (dM), F(X) S R. 

Similarly we decompose fj c (X) into a tangential part G(X) and a normal part F{X). 

Finally, let S : T(dM) -> T(dM) denote the shape operator defined by S(X) = 
V^n where n is the inward unit normal to dM. Then S is a self-adjoint opera- 
tor whose eigenvalues are the principal curvatures of dM, trS — H is the mean 
curvature of dM and det(S) = K ext is the extrinsic curvature of dM. (Note that 
with our sign convention for S 1 , the principal curvatures are positive when dM is 
concave.) 

With this notation established, we make the following computation: 

Lemma 3.1. Let a, a be the 1-forms on dM dual to the vector fields v, v on dM. 
Then, using the notation defined above, G,G : T p (dM) — » T p (dM) and F,F : 
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T p (dM) — ► R are given by 

(12) G = Dv + hS -hJ, G = Dv + hS + hJ 
and 

(13) F = dh- Sa - *or, F = dh- Sa + *a 

where D is the exterior covariant derivative on dM , S is defined on 1- forms a by 
Sa(X) = Sv ■ X = a(SX), and J : T(dM) -> T(dM) is the rotation by tt/2 given 
by JX = n x X (so J 2 — —identity). 

Proof. For X £ T(dM) we have, using equations © and l[TU]). 

T] C (X) — Vxu — uxX 

= Vx(v + hn) - (v + hn) x X 

= V x v + hVxn + X(h)n-v x X -hnx X, 

Thus the tangential part G : T p (dM) -> T p (dM) is given by 

G(X) = V x v + hS(X) -hnx X = (Dv + hS - hJ)(X) 

where V, D denote the Riemannian connection and exterior covariant derivative on 
dM. Further the normal component F : T(dM) — > R is given by 

F(X) = r) c (X) ■ n = X(h) -v x X -n+ {V x v) -n 

= grad h ■ X — n x v ■ X — v ■ Vjn 

= (gradh-Jv-S(v))-X 

= (dh-*a- Sa)(X). 

Similarly, we find fj c = Du + rot u has tangential part G : T p (dM) —* T p (dM) 
given by 

G{X) = V x v + hS(X) +hnxX= (Dv + hS + hJ){X) 
and normal part F : T(dM) — > R given by 

F(X) = X{h) +vx X-n+ (Vx«) • n = (dh + *<r - Sa)(X). 

□ 

We can view F, F and G, G as real- valued and vector-valued 1-forms on dM. It 
is then possible to define the wedge products F A F and GAG, where in the latter 
case we also use the dot product on TM from the hyperbolic metric on M to obtain 
a real-valued 2-form. Then our boundary term can be expressed in terms of these 
wedge products as 

b = / rjc A fj c = / FAF+GAG. 

JdM JdM 

We will now compute the two summands in this expression separately. First re- 
call that if to is a 1-form on a Riemannian 3-manifold M, then its exterior derivative 
duj satisfies 

(14) du(X,Y) = Xuj(Y) -Yu(X) - w([X,Y]) 

for all vector fields X, Y on M. We can also rewrite this using covariant derivatives: 

(15) duj(X,Y) =V x lu(Y) -Vy^(X) 
since V X Y - V Y X = [X, Y). 
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To analyze the boundary term J dM GAG we will use the following. 

Lemma 3.2. Let v,v be vector fields on dM with dual 1- forms a, a, and let dA 
denote the area 2-form on dM . Denote by 6 the adjoint of the exterior derivative 
on dM . Then 

(1) J A S = tr S dA = H dA where H is the mean curvature of dM, 

(2) J A Dv = — * 8a, 

(3) Dv A S = d(Sa), 

(4) S A S = 0, 

(5) J A J = 0, _ 

(6) J gM Dv A Dv = (Ka,*a), where K is the Gaussian curvature of dM and 
(•, •) is the L 2 inner product on dM. 

Proof. Let e±, e^ be an oriented orthonormal basis for T p (dM). Then for any linear 
operator L : T p (dM) -> T p (dM) we have 

(JAL)(ei,e 2 ) = J(ei) • L[e 2 ) - J(e 2 ) • £(e x ) 

= e 2 ■ L(e 2 ) + ei • L{e{) = tr L. 

Hence J A 5 = tr S dA = H dA and J A £>v = tr Dv dA = div v dA = — * 5a. 
Next we prove part (3). For X,Y E T p (dM) we have 

(DvAX)(X,Y) = V x v ■ S(Y) - V Y v ■ S(X) 
= II(V x v,Y)-II(V Y v,X) 
= XII(v,Y)-II{v ) VxY)-{V x n){v ) Y) 

- (YII(v, X) - II(v, Vyl) - (V Y H){v, X)) 

where II(X,Y) = S(X) Y = X ■ S(Y) is the second fundamental form. But 
(\7 xII)(v,Y) = (VyII)(v,X) by the Codazzi-Mainardi equations for a space of 
constant curvature (see, for example, p2, chap.l, Thm 11 and Cor 12]). Hence 

(Dv A S)(X, Y) = XII(v,Y)-YII(v,X)-II(v,[X,Y}) 

= X(Sv ■ Y) - Y(Sv -X)-Sv- [X, Y] 

= X(Sa(Y))-Y(Sa(X))- Sa([X,Y}) 

= d(Sa)(X,Y). 

Finally, we have S A S = and J A J = by the skew-symmetry of the wedge 
product, and for vector fields X, Y on dM we have 

D 2 v(X,Y) = ViVy« - VyVxll - V[ X ,Y]V = R{X,Y)v, 

where R is the Riemann curvature tensor on dM. But R(e\ 1 e2) is infinitesimal 
rotation by —K, hence D 2 v = —KJv dA where dA is the area 2-form on dM and 

/ Dv A Dv = / d(y A Dv) - vA D 2 v = Kv ■ Jv dA = (Ka, *a) 

JdM JdM JdM 

by Stokes' theorem. This completes the proof of the Lemma. □ 

Using this result and Lemma |3. II we obtain 
Lemma 3.3. 

/ GAG = (h,*dSa)-{h,*dSa)+{h,Sa)+(h,6a)-({Hh,h)+(Hh,h))+{Ka,*a), 

JdM 

where (a, (3) = J gM a A *(3 denotes the L 2 inner product on dM. 



THE SHAPE OF HYPERBOLIC DEHN SURGERY SPACE 



15 



Next we study the boundary term /„„ F A F = —{F, *F). From Lemma l3~Tl we 
have 

F = dh — So — *o and F = dh — So + *o. 
Using this we obtain 

/ F hF = -(dh - So - *o,*dh- *So - o) 

JdM 

= - (dh, *dh) + (dh, *So + o) + (So + *o, *dh) - (So + *o, *So + o). 
But the first term vanishes since 



-(dh,*dh) = / dhAdh = d(h A dh) 

JdM JdM 



= 0, 



so we obtain 
Lemma 3.4. 



/ F A F = (dh, *So + o) + (So + *dh) - (So + *a, *So + o) . 

JdM 

Combining the previous results gives the following 

Theorem 3.5. Let s = u + iu be a section of the bundle E and let dES — r\ c + ifjc 
denote its image under the coboundary operator dgs. Suppose that u — v + hn, 
u = v + hn are the decompositions into tangential and normal parts of the vector 
fields u,u. Denote by o,o the l-forms on that are dual on dM to the vector fields 
v,v, respectively. Then the Weitzenbock boundary term equals 

n c Afj c = 2((h,5o) + (h,So)) - (So + *o,*So + o) 

dM 

-((Hh,h) + (Hh,h)) + (Ko,*o). 
Proof. From the previous lemmas we have 



n c Afj c = / FAF + GAG 

dM JdM 

= (dh, *So) + (dh, o) + (So, *dh) + ( * o, *dh) 

— (So + *cr, *So + o) 

-(h, *dSo) + (h, So) + (h, *dSo) + (h, 8a) 

- ( (Hh, h) + (Hh, h) ) + (Ko, *a) . 

We can simplify this sum by noting that (dh,o) — (h,So), ( * o,*dh) — (h,5o), 
(dh, *So) — (h, *dSo) — 0, and (So, *dh) + (h, *dSo) = 0. This gives the result. □ 



The computation of the boundary term in Theorem 13.51 is valid for general el- 
ements rj c ,fj c £ 17 1 (M;TM) that are of the form rigs = n c + f\ c i. However, this 
boundary term is primarily of interest when it comes from the Weitzenbock formula 
J7]), as discussed in the previous section. Specifically, we are interested in the case 
when, in a neighborhood of the boundary, both r\ c ,f] c are symmetric and traceless 
when viewed as elements of Horn (TM,TM). Writing s = u + iu as before, the 
condition that n c and f) c are symmetric is equivalent to the equations (derived in 
QUI Section 2]): 

(16) Iu — curl u, 2u — — curlu 
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Here we view the curl of a vector field in 3-dimcnsions as itself being a vector 
field. Then, on the boundary of M, the normal component of the 3-dimensional curl 
of u is just the (scalar) 2-dimensional curl of v, the tangential part of u. A similar 
statement holds for the normal component of curlu. Since the normal components 
of u, u equal h, h, respectively, we obtain: 

(17) 2h = curlu = *da and 2h = — curlw = — * da, 

where a, a are the 1-forms on dM dual to v, v, respectively, and d denotes the 
exterior derivative operating on forms on dM . 

Note that the equations (|16p only hold in a neighborhood of the boundary. 
However, since all our computations are local to the boundary, this will suffice. It 
turns out that (fl7J|) also implies that i] c and fj c are traceless (using divcurl = 0), 
but we will not use this in our computations. 

Using (|17j) we can rewrite b in terms of a and a. We compute 

2 (ft., 5a) = — ( * da, da) = — (8(*a), 5a) = — ( * a, d5a) , 
2(h, 5a) = (* da, da) = — (da, d(*a)) = — (Sda, *a) , 

and note that 

( * da, *da) = (5(*a), 5(*a)), ( * da, *da) = (da, da). 
Then we obtain 

(18) b = -(* a,d6a + 5da) - (Sa + *a,*Sa + a) 

--((H6(*a),8(*a)) + (Hda,da)) + (Ka,*a). 

Equations (|16[) provide relations between u and u, equations (|17p coming from 
the normal component of those relations. Similarly, the tangential component of 
(|16[) implies that v, hence a, can be expressed in terms of u and its derivatives 
in a neighborhood of the boundary. In the next section wc will define boundary 
conditions that will allow us to express a on dM purely in terms of a and its 
tangential derivatives. In particular, the expression will not involve the normal 
component h or normal derivatives of a. 

There are many possible boundary conditions of this sort and we denote by A 
a general linear differential operator on 1-forms on dM. Then, as a simplifying 
notational device, we can express the relation between a and a on dM as: 

(19) a = *Aa 

Finally, it is useful to define a linear operator S on 1-forms by S — — * 5*, so 
that *Sa = S * a. It is easy to check that this operator satisfies 

SS = det(S)J = K ext I, 

where K ext is the extrinsic curvature of dM. It follows, since the curvature of the 
ambient hyperbolic manifold equals —1, that 

§S-I= (K ext - 1)1 = K mt I = KI. 

As before K denotes the intrinsic Gaussian curvature of dM. 
Rewriting equation (|18p using the operators S and A we find 

6 = (Aa,Aa) + ((A ~ S)a, (I - SA)a) 

--((H5Aa,6Aa) + (Hda,da)) - (Ka,Aa). 
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Hence we obtain 

Theorem 3.6. If 2h = curlu = *da, 2h = — curl 6 = — * da, S — — * S* and 

a = *Aa, then the Weitzenbock boundary term is given by 

(20) b = ((A-S)a ) (I-SA)a) + (Aa ) (A-K)a)-~((H5Aa,5Aa)+(Hda > da)). 

We observe that the expression for the boundary term in Theorem 13.61 can be 
viewed as a quadratic form on 1-forms a. Except for the operator A, the basic terms 
in this quadratic form come from the geometry of the boundary of M. In particular, 
S, K, H are the shape operator, Gaussian curvature and mean curvature of dM. 
So, given the manifold M with its boundary dM, the only flexibility we have on 
this quadratic form is the tangential operator A. We can attempt to control this 
operator by our choice of boundary conditions when solving equation ([3]). Recall 
that the 1-form r in ([3]) is dual to the vector field u which in turn determines u by 

Our goal is to find boundary conditions which determine an operator A with the 
property that this quadratic form is positive semi-definite; i.e., so that the boundary 
term (l2Cfl) is non-negative for all a. 

We now specialize to the case when M has tubular boundary. Then each bound- 
ary component is topologically a torus, and the principal curvatures fci,fc 2 are 
constant so that the mean curvature H = k\ + fc 2 is constant. The extrinsic curva- 
ture is K ext — k\k2 = 1 and the intrinsic curvature is K — 0; i.e., the torus is flat. 
The operator S equals 5 1-1 in this case and both S and S are parallel. Then the 
boundary term simplifies to 

(21) b = ((A - S)a, (I - SA)a) + (Aa, Act) 

- — ({SAa,SAa) + (da, da)). 

If we denote by Aq the 0-th order part of the operator A (i.e., the part that 
involves taking no derivatives), then the 0-th order part of this quadratic form is 
simply 

((A Q - S)a, (I - SA Q )a) = -{(A - S)a,S- 1 (A Q - S)a), 

where we have used that fact that S — S^ 1 to obtain the second expression. Since 
S 1-1 is a positive operator this quantity is non-positive. Our only hope of having a 
non- negative quadratic form is to choose A so that Ao — S. 

The computations below show that the choice of A — S in fact does lead to a 
non-negative quadratic form. However, as will be discussed in the next section, we 
have been unable to find an elliptic boundary value problem that leads to this value 
of A. Nevertheless, we are able to find such a boundary value problem that leads 
to a slightly perturbed value of A that still defines a non-negative quadratic form. 

Suppose the tangential operator A equals 

(22) A = S+^6d, 

where e > is a constant. Thus, in the above Weitzenbock boundary term we have 
A - S = §<M and / - SA = I - SS - %S8d = -%S6d. 

In the next section we will show that it is always possible to solve equation ([3]) in 
such a way that the 1-forms a, a arising from the solution satisfy the relation (Q3j|) 
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with this value of A. For now, we will assume that this can be done and complete 
the computation of the boundary term with this value of A. 

Now 5 A = 5S + ^SSd = SS and (5da,d5a) = (da,ddSa) = 0, so equation (12"T]) 
becomes 

T] c Afj c = (Sa, Act) - — ((5Sa,SSa) + (da, da)) 

dM 4 



(I) (Sda, SSda) + | (5da, Sda) . 



and we want to find geometric conditions on dM guaranteeing that this boundary 
term is non-negative. 

First we consider (SSa, SSa) + (da, da) and (Sa, Aa) . Since the metric on dM is 
Euclidean, we can choose a parallel orthonormal frame field e\, e 2 on dM consisting 
of eigenvectors of S with eigenvalues k\,k 2 at every point. Let 9\,0 2 be the dual 
1-forms on dM and write a = a\B\ + 02 #2- Now, using equation (fTS)) . 

da = (V1CT2 — V2CT1) dA and - 5(Sa) = (k x Vidi + fc 2 V 2 cr2) dA 

where V, = V e; for i — 1,2. Hence 

(*SSa) 2 + {*da) 2 = ^(VktO 2 + k 2 {X7 2 a 2 ) 2 + {X7 1 a 2 ) 2 + (V 2 cti) 2 

+2 (fcife2(Vio-i)(V 2 (r 2 ) - (Vi<r 2 )(V20-i)) . 

Since k\k 2 = 1 the bracketed terms become 

V\aiV 2 a 2 ~ V\a 2 V 2 ai = da\ A da 2 (e\,e 2 ), 

and their integral over dM is 



/ dai A da 2 = / d(ai A da 2 ) = 0. 

JdM JdM 



Hence 

\\SSaf + \\da\\ 2 = ^HVkhH 2 + k 2 \\V 2 a 2 \\ 2 + ||V 1( i2|| 2 + ||V 2 a 1 || 2 
where || • || is the L 2 -norm on dM. Using integration by parts, 
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(Sa,Aa) = / Vfc^Acr; = V fcj gradaj 2 = V h\\ V^if. 

Combining the last two equations and using kik 2 = 1 we obtain 
A(Sa,Aa) - H(\\SSa\\ 2 + \\da\\ 2 ) 

= (3 - fc^MVifTiH 2 + (3*1 - fc2)||V 2f Ti|| 2 + (3fc 2 - fci)||Via 2 || 2 + (3 - k 2 )k 2 \\S7 2 a 2 \ 
2 

= (3-fc 2 )^||V^-|| 2 

To examine the other terms, write Sda = ai#i + a 2 9 2 . Then 

SSda = k 2 ai9i + k\a 2 9 2 
since S = S^ 1 and k\k 2 = 1. Hence 

^(Sda, Sda) - (|) (Sda, SSda) = ~J ((a 2 + a 2 ) - |(fc 2 a 2 + half) dA 

((1 - £ -k 2 )a\ + (1 - £ -h)al) dA. 



DM 
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This will be non-negative provided — < — = k\ and — < — = fej that is, if 

2 K2 2 fci 

< — < min(fci, 1*2) ■ 

This gives our final conclusion: 

Theorem 3.7. Let M be a hyperbolic 3-manifold with tubular boundary and let 
A = S+^5d where e > is a constant. If2h = curl?; = *do~, 2h = — curlw = —*da 
and a = *Ao~, then the Weitzenbock correction term b — J dM r\ c A fj c is 
2 

- ^(S-^HV^-f + l / ((^2 ~ + (fci - |)al) dA. 
i,j=i -' 9M 

Hence the boundary term is non-negative if the principal curvatures fci,fc2 satisfy 

— < ki < k 2 < V3 

and e < 2k\ . 

4. Boundary Values 

In this section we will describe a boundary value problem that will allow us to find 
harmonic representatives for infinitesimal deformations of hyperbolic 3-manifolds 
with tubular boundary whose boundary values are of the form discussed in the 
previous section. This will allow us to make statements about the boundary term 
in the Weitzenbock formula which, in turn, will lead to local rigidity results for such 
3-manifolds. Those results and other applications will be discussed in Section 5. 

In Section 2 we saw that finding a harmonic representative uj for an infinitesimal 
deformation amounts to finding a real-valued 1-form r which is a solution to the 
equation (A + 4)t = (. Here, £ is a smooth, real-valued 1-form which equals zero 
in a neighborhood of the boundary and A is the usual Laplacian on real-valued 
1-forms on M. The 1-form r is dual to a vector field u on M which is the real 
part of an E- valued section s and the coboundary c?^ s is added to the original E- 
valued 1-form in order to make it globally harmonic. The boundary behavior of r 
determines that of s and hence of g?£ s, providing information about the boundary 
values of u). 

In order to have any control over the behavior of r near the boundary, it is 
necessary to put restrictions on the domain of the operator (A + 4). However, the 
restrictions must still allow the above problem to be solvable. Below, we will define 
boundary data that the real- valued 1-form r must satisfy which make this operator 
elliptic, self-adjoint with trivial kernel. Standard theory (Chapter X in [T3]) then 
implies that the above problem is uniquely solvable; when £ is smooth, as it is in 
our situation, the solution t will be smooth. 

There are many choices for such boundary conditions. Standard examples in- 
clude prescribing that either the value or the normal derivative of r be zero, anal- 
ogous to Dirichlet and Neumann conditions for the Laplacian on real-valued func- 
tions. However, our choice is motivated by the further condition that the resulting 
Weitzenbock correction term b defined in (fTTj) be positive. None of the more stan- 
dard choices of boundary data have this property. 

In order to describe our boundary conditions we first need to establish some 
notation. 
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The above Laplacian on 1-forms, A, equals 

dS + Sd 

where d is exterior differentiation on M and 5 its adjoint. We will denote by d and 
(5 the corresponding operators on dM . Similarly, we use the notation * to denote 
the 3-dimensional Hodge star operator on forms, reserving the notation * for the 
corresponding operator on the boundary. 

It is also useful to define operators ds = Sd and Ss = 8S operating on func- 
tions on the boundary and on tangential 1-forms respectively. Here S is the 2nd 
fundamental form or shape operator on the boundary, with normal chosen so that, 
in our situation with concave boundary, S is positive definite. We also define 
As = Sd + dsSs, which acts on tangential 1-forms. 

Recall the basic setup from Section 2: We begin with an E- valued 1-form 
which represents the cohomology class in H 1 (M ; E) determined by our infinitesimal 
deformation. It satisfies d^G) = 0, but, in general, 5e& ^ 0. To find a harmonic 
representative we must find a globally defined ©-valued section s satisfying the 
equation SecLes = — #e<2>. Then u — Q + rigs is a harmonic representative in the 
same cohomology class as Cj. 

Decomposing s into its real and imaginary parts, we write s = u + i u where u 
and u can be viewed as vector fields on M. We can assume that 2u = curlu (by 
choosing s to be a canonical lift, see Section 2 of [TO]. As discussed in Section 3 
of the current paper, this is equivalent to the real part of rigs being symmetric.) 
Thus, it suffices to find u. This is equivalent to solving the equation (A + 4)r = £ 
where £ is the 1-form dual to the real part of the ©-valued section — 5e&> and r 
is the 1-form dual to u. The equation 2u = curlu is equivalent to the equation 
2f = Sdr, where f is the 1-form dual to u. 

As we have done before, we can decompose 1-forms on (a neighborhood of) the 
boundary into their normal and tangential parts. In particular we write 

(23) t = hdr + a and 2f = *dr = 2(h dr + a), 

where dr denotes the 1-form dual to the inward pointing unit normal and <r, a are 
tangential 1-forms. 

We now describe a 1-parameter family of boundary conditions, parametrized by 
a parameter e. It is assumed that e > and is a constant. Using the notation 
established above, the boundary conditions can be expressed as: 

(24) It - 2e{S s a - 2h) = 

(25) 2a-*(2Sa + e(A s a-2d s h)) = 0. 

At the end of this section we will show that the boundary value problem of 
solving (A + 4)t = C subject to these boundary conditions is elliptic and that the 
operator is positive, self-adjoint. This implies that there will be a unique solution 
and that the solution will be smooth on the entire manifold with boundary. 

However, in order to provide some motivation for choosing these fairly compli- 
cated boundary conditions, we will first assume the existence of such a solution and 
analyze the properties of the harmonic -©-valued 1-form to that we obtain from t. 

Note that w is dE closed and, hence, is the image under dE of a locally defined 
section of E whose real part is a locally defined vector field. The divergence of 
this vector field is just the trace of the real part of u>, viewed as an element of 
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Horn (TM,TM). Thus, although the vector field is only locally defined, its diver- 
gence is a globally well-defined function. By abuse of language we will refer to this 
as the "divergence of u>" . 

The main step is to show that when r satisfies the above boundary conditions, 
the resulting harmonic _E-valued 1-form 10 has divergence identically zero. This will 
imply that the stronger harmonicity equations (j4]) and hold and, thus, that the 
results from Section 2 and the computations from Section 3 all apply. It will also 
show that t in fact satisfies boundary conditions that are stronger and simpler than 
HMD and (USD. 

Proposition 4.1. Let lo = lo + cIes be a harmonic E-valued 1-form on a compact 
hyperbolic 3-manifold with tubular boundary, where Co is in standard form near the 
boundary and s~u + iuisa global section of E. Let t,t be the 1-forms dual to 
the vector fields u, u, respectively. If r and f are decomposed as in (23\) and satisfy 
the boundary conditions \21$ and i25\) . then the divergence of lo is identically zero. 

Proof. We denote by tr the divergence of the harmonic deformation lo. Harmonicity 
of lo implies (see [10l Lemma 2.4)] that 

(<ta+4) tr = 0. 

This equation holds on all of M and, taking the L 2 dot product on M of tr with 
this equation, we conclude that ((8d + 4) tr, tr) = 0. Integrating by parts gives 

(d tr, dtr) + 4(tr, tr) + / trAidtr = 

JdM 

where the boundary is oriented using the inward normal. If we show that the 
boundary integral 

/ trAidtr 
JdM 

is non-negative, it will follow that tr — and dtr = on all of M. In particular, 
we will have shown that lo is divergence-free. 

Since Cj equals some standard harmonic _E-valued 1-form in a neighborhood of 
the boundary, and since standard forms are all divergence-free, the divergence of 
lo = Cj + dss just equals the divergence of rigs in a neighborhood of the boundary. 
By definition the latter equals the divergence of the vector field u which is the real 
part of s; this equals —St, since r is the 1-form dual to u. 

Since tr = — St in a neighborhood of the boundary, we can use the boundary con- 
dition (|2"4"|) on t when computing the boundary integral. Also, since Cj is harmonic 
in a neighborhood of the boundary, 8eCj = in a neighborhood of the boundary 
which in turn implies that SEdgs is zero near the boundary. As discussed in Section 
2, this means that (A + 4)r = near the boundary. 

The second term in the integrand becomes — *cL6t, but, because r satisfies (dS + 
8d + 4)r = 0, this equals *(5d + 4)r. Since the integral is over the boundary, only 
the tangential part of the integrand appears. The tangential part of *5dr = d*dT 
equals d(2a) where 2tr is the tangential part of idr and d is the exterior derivative 
on the boundary. The tangential part of *r equals h dA, where dA is the area form 
on the boundary. Rewriting 2da as 2S(*a) dA, the boundary integral can then be 
written as 

[ (-5r)(4/i + 8(*2a)) dA. 

JdM 
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Using |24|) and (|25|). the boundary integral becomes 

/ 2e(2h - 8 s a){Ah - 8{2Sa + e(A s a - 2d s h)) dA. 

JdM 

This simplifies to 

2e I (2h-S s a) [{2 + eSd s )(2h - S s a)) dA, 

JdM 

which equals 

2e ((2/i - Ssa), (2 + eSd s ){2h - 5 s a)) aM , 
where (•, -)qm denotes the L 2 dot product of 1-forms on dM. 

We see that 2 + eSds = 2 + eSS^S^d is a positive operator since S is positive 
self-adjoint and hence has a positive self-adjoint square root. Thus the boundary 
integral is non-negative and we conclude that the deformation is divergence free on 
M as desired. □ 

The fact that u> is divergence-free provides us with further information about 
the boundary behavior of the solution r. In particular, since tr = — St near the 
boundary, one concludes, as in the previous proof, that St = near the boundary. 
Together with the first boundary condition (|2"4")l . this implies that 

2h = 5 s a, 

and the second boundary condition (f2"5|) simplifies to 

2a = *(2Sa + eSda). 

The computations of the general Weitzenbock correction term in the previous 
section were all expressed in terms of a general differential operator A on tangential 
1-forms. The operator A expresses the relation between the two tangential 1-forms 
a, a and is determined by the equation (fl9|) a — *A<r. With our choice of boundary 
conditions ([24]) and (|25[) . we see that 

Act = Sa + —eSda 

which is precisely the value (|2"2"|) for the operator A that we wished to obtain. It 
was with this result in mind that we were led to our boundary conditions. 

We noted in the previous section that in order for the Weitzenbock correction 
term b in (|2"0|) to be positive, it is necessary for the 0-order term of Aa to equal 
Sa and that choosing Aa — Sa did give a positive value for b. Indeed, a natural 
choice for our boundary conditions would have been to set e = in (|2"4")) and ([23)1 : 
this would have led to the value Aa = Sa. However, those conditions do not lead 
to an elliptic boundary value problem and it was necessary to perturb this natural 
choice to obtain an elliptic problem. To do so in such a way that the resulting 
harmonic deformation was divergence-free and so that b was still positive required 
some delicacy and led to the more complicated form of the boundary conditions. 

We are now in a position to apply the conclusions from the previous sections. 

Let uj = uj + dES equal the harmonic E- valued 1-form obtained from our bound- 
ary value problem, and decompose u into its real and imaginary parts as r\ + if} 
where r\, r) are elements of Horn (TM, TM). We have assumed that r\ is symmetric 
and Theorem 14.11 implies that it is traceless. Therefore the stronger harmonicity 
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equation j4| holds. As discussed in Section 2, this implies that fj is also traceless and 
symmetric and equals *Drj. This, in turn, allows to conclude that the Weitzenbock 
formula ([7]) holds. 

In a neighborhood of the boundary, we can write w = ojq + lu c , where ojq is a 
standard form and oj c = d^s. Since the real and imaginary parts of ojq are also 
symmetric and traceless, the same will be true for oj c . Writing uj c = d^s = r) c + ifj c , 
we see that the hypotheses of Theorem 13.61 hold for r) c ,fj c . Since we have also 
concluded that, for any oj obtained using our boundary conditions, the operator A 
satisfies (|22[) , Theorem 13 . 71 applies . If we write ojq = r)o + ifjo, then, as discussed in 
Section 2, fj = *Drj and fjo = *Drjo. It follows that fj c — *Dr] c as well. Hence, the 
boundary integral b in equals —b R (r) c ,r] c ) where b R (r) Cl r) c ) is the contribution 
from the correction term oj c to the Weitzenbock formula ([7|): 

b= I r)c^V c = T] c /\*Dr) c = -b R (ri c ,ri c ) 

JdM JdM 

Thus, assuming the existence of the solution r of our boundary value problem 
(proved in Theorem 14.41 below) . Theorem 13 . 71 implies : 

Theorem 4.2. Let M be a hyperbolic 3-manifold with tubular boundary whose 
principal curvatures k\,k 2 satisfy ^= < k\ < ft2 < Then, for any infinitesi- 
mal deformation, there is a harmonic representative oj so that the correction term 
bu(j] c ,T] c ) is non-positive. 

Recall from Section 2 that, for any divergence- free harmonic oj, we have the 
Weitzenbock formula (0 [8]) : 

I \Dv\ I 2 + I hi I 2 = b R (ri,Ti) =b R (r ]0 ,T l0 ) + b R (r )c ,r} c ). 

We immediately obtain the following corollary which will be crucial in proving the 
applications in the next section. Note that the standard form ojq depends only 
on the infinitesimal variation of the holonomy of the boundary; thus, it and its 
real part rjQ are invariants of the cohomology class of the infinitesimal deformation, 
independent of the choice of representative. 

Corollary 4.3. Let M be a hyperbolic 3-manifold with tubular boundary whose 
principal curvatures fci,fc2 satisfy ^= < k\ < k 2 < \A3- Then, for any non-trivial 
infinitesimal deformation of M we obtain b R (r]o,r]o) > 0. 

Finally, we must justify our claim that we can always solve our given boundary 
value problem. 

For our purposes a differential operator on a manifold M with boundary consists 
of a differential operator P from C°°(F), the C°° sections of a bundle F over M 
to C°°(G), where G is another such bundle, together with a collection {pi, ...,p r } 
of differential operators from C°°(F) to (&iC°°(Gi) where (BiGi is a direct sum 
decomposition of the bundle G, restricted to the boundary. Such an operator will 
be denoted by (P; {pi, ...,p r }). In our case, the bundles F and G are both equal to 
the bundle of 1-forms on M and we decompose this bundle on the boundary as the 
direct sum of its normal and tangential parts. The main operator P is A + 4 and 
Pi,P2 equal the operators on the left-hand sides of (f2~4|) and (j25|) . respectively. 

The remainder of this section will be devoted to proving the following theorem. 
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Theorem 4.4. Let M be a compact hyperbolic 3-manifold with tubular boundary. 
Then, for any constant e > 0, the differential operator (P;{pi,P2}) on the bundle 
of C°° real-valued 1- forms on M defined, using the notation above, by 

P{r) = (A + 4)(r) 

pi(r) = 5(r) - 2e(S s a - 2h) 

P2 (t) = 2a-*(2Sa + e(A s a-2d s h)) 

is elliptic. On the subspace where pi(r) = P%(t) = it is positive and self- 
adjoint. In particular, for any smooth 1-form there is a unique solution to 
P(t) — CiPi( T ) — P2(t) — and that solution is smooth on all of M . 

Proof. We will first show that this operator is self-adjoint and positive. 

Recall that A = d5 + 5d. Then for any real- valued 1-forms T,ip on M, integration 
by parts gives us: 

(At, tp) = {dr, dyj) + {St, 5yj) + B(t, i/>), 

where (■, ■) is the L 2 inner product on M and B(t, ip) is a boundary term which is 
given by an integral over the boundary. 

The operator A (hence A + 4) is self-adjoint as long as 

(Ar,^) - (A^,r) = B{r^) - B(J>,t) = 0. 

The operator A + 4 will have trivial kernel as long as 

((A + 4)t,t) >0 

for any non-zero r. Letting r = ip above this will be guaranteed as long as we have 

B(t,t) > 0. 

Using Green's identity, we obtain the following formula for the boundary term, 
where the boundary is oriented with respect to the inward normal: 

B(T,ip) = — ( / *dr A V' + 5rA*ip). 

JdM 

Again we use the notation * to denote the 3-dimensional star operator on forms, 
reserving the notation * for the corresponding operator on the boundary. 

As before we decompose r as t = hdr + a and let 2a equal the tangential part 
of *dr. If we decompose ip as ip — k dr + (f>, we can write 

B(r,ip) = / 2<j) A a - St A *k. 

JdM 

Using the boundary conditions (f24|) and (|25|) . the boundary term becomes 
B{t,i/j) = I (/>A*(2Sa + s(A s a-2d s h) + 2e{2h - S s cr) A *k) 

JdM 

= 2(4>, Scr) + e{{4>, Sdo- + d s (5 s c7 - 2h)) + (-2*, {S s a - 2h))) 

= 2(0, Sa) + e((#, da) + ((6 s <t> ~ 2k), (5 s a - 2h))) 

where (•, •) denotes the L 2 inner product on forms on the boundary oriented by the 
inward normal and we have used the definition As = Sd + dsSs- 

It is apparent from this formula and the fact that S is symmetric and positive 
definite, that the boundary term is symmetric in r and i\) and non-negative when 
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i/j = t. It follows that the operator A + 4 is self-adjoint and positive definite with 
these boundary conditions. 

Finally, we must show that this boundary value problem is elliptic. To see that 
the boundary conditions lead to an elliptic boundary value problem, we consider 
each of the operators, P, pi,P2 where the ranges of the two latter operators are the 
sub-bundles of normal and tangential parts of 1-forms on the boundary. We then 
take the top order terms of each of these operators. It is a subtlety of differential 
operators on bundles that ellipticity may depend on the choice of decomposition of 
the target bundle, since this affects what the top order terms are. 

To show that the system is elliptic one considers the symbols of the operators. 
This amounts to looking at the system in local coordinates, fixing the coefficients 
of the operators by evaluating at a boundary point, and taking only the top order 
terms in each operator. One then considers the homogeneous, constant coefficient 
problem in the upper half space of 1" given by these simplified operators. We refer 
the reader to [13], Chapter X or [2], Appendix I for a full description. Below we 
will see how the process works in our specific case. 

For A + 4 taking the symbol simply gives the standard laplacian in M 3 which is 
well-known to be elliptic. Since 2a equals the tangential part of *dr, it is obtained 
by applying a first order operator to r. One then easily sees that pi is of first order 
and P2 is of order 2. Taking the top order terms, the boundary operators simplify 
to 

St - 2e(5 s a), 
eA s cr. 

These operators are still defined in terms of the hyperbolic metric. It is easy to 
check that, taking natural orthonormal coordinates at any point on the boundary 
torus the coefficients of the top order terms are independent of the point chosen. 
The operators simply become the same operators viewed in the upper half-space 
K.+ = {(x,y,t)\t > 0} with the standard Euclidean metric. Multiplication by S 
becomes multiplication in the (x, y)-plane by the diagonal matrix with diagonal 
entries k = ki, k^ 1 — ft2 which are the principal curvatures of the boundary torus. 
We are now left with the simplified system of solving At = in the upper half-space 
with homogeneous boundary conditions determined by these simplified boundary 
operators. By definition, the original system is elliptic if and only if this simplified 
system has no non-trivial bounded solutions. It suffices to show that there are no 
non-trivial bounded solutions r(x, y, t) of the form 

where £ = (a, b) is any non-zero vector in the boundary plane. The solutions, 
r(i), to At = for a given choice of £ are linear combinations of e^'* and e - ' 1 "'*. 
Since we are only interested in bounded solutions, only scalar multiples of the latter 
function will appear. In particular, we have that r'(t) = — |£|r(i). 

We decompose r into its normal and tangential components which we again 
denote by h and a, respectively. Viewing r as a 3-dimensional vector field and a as 
a 2-dimensional one, —8t is the divergence of r which equals hi plus the divergence 
of a. Similarly, — Ssct is the divergence of Sa. Note that all calculations are done 
with respect to the Euclidean metric. 
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For solutions of this form with £ fixed, the boundary conditions become 



fe(0)|C|-iC-*(0) 



i2s(-Sa(0) = 
eLs(()a(0) = 0, 



where Ls(C) is the matrix computed below. 

We'll see that Ls(Q is invertible, so the second boundary condition implies that 
<r(0) = 0, and thus, the first implies that h(0) = 0. This means that r(0) = and, 
since r'(0) = — |£|r(0), we conclude that any solution must be trivial. 

To compute the matrix in the second boundary condition, we view a as a 1-form 
on the boundary of the upper half-space. Recall that Ag = 5d + dsSs = Sd+ SdSS . 
It's a standard calculation that for £ = (a,b), the symbols for Sd and dS are, 
respectively, multiplication by the matrices 



—ab 



-ab 



and 



a 2 ab 
ab b 2 



The matrix S is diagonal with entries k, k 1 on the diagonal. It follows that the 
matrix corresponding to the operator is 



MO = 



k 2 a 2 









- At 



2 h 2 



which is clearly invertible for all (a, b) ^ 0. 



□ 



It is worth pointing out that, when e = in the original boundary conditions, 
the top order term in the second equation is 2<r which is of first order and the 
simplified conditions become 

St = 
2ct = 0. 

An easy calculation shows that the first order part of 2tr equals a' — dh, so that 
the corresponding linear equations used to determine ellipticity are 

fc(0)|fl-iC-*(0) = o 
a(0)|C|+i/i(0)C = 0. 

This system is seen to have a non-trivial solution given by icr(0) = /i(0)t|t for 
any £ ^ 0, and thus the system is not elliptic. This is the reason we needed to 
perturb the system by adding a small second order term in the second boundary 
condition. The term added to the first boundary condition was necessary to keep 
the system self adjoint. The precise form of these added terms was determined by 
other conditions necessary to conclude that the system was divergence free and led 
to a positive value for the Weitzenbock boundary correction b. 



5. Applications to hyperbolic Dehn Surgery Space 



We now apply our harmonic deformation theory to study generalized hyperbolic 
Dehn surgery as introduced by Thurston in [18] . 
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5.1. Geometry of tubular boundaries. 

Let M be a compact orientable hyperbolic 3-manifold M with tubular boundary 
and let T = Tr be a torus in dM of tube radius R < oo. Then T has principal 
curvatures k\ = cothi? and &2 = tanhi?, and the intrinsic Euclidean metric on T 
has the form E 2 /r where r = Z 2 acts as translations of E 2 , and the holonomy of 
the Euclidean structure gives an isomorphism h : tt±(T) — > T. 

The developing map for M restricts to an isometric immersion $ : T = E 2 — > H 3 
taking the universal cover T of T to the surface of a cylinder of radius R in H 3 ; <!> 
is uniquely defined up to composition with isometries of H 3 . 

Explicitly, if we choose standard Cartesian coordinates on E 2 with x\,X2 coor- 
dinate axes in the directions of the principal curvatures k\ , k% respectively we can 

take <fr(xi,X2) to be the point with hyperbolic cylindrical coordinates 9 = — : — 

sinn it 

and C = • Then for each 7 g -k\ (T) , the complex length of 7 is given by 

coshi? 

(26) £(7) ''- ■ " 



cosh R sinh R 

where the Euclidean holonomy ^(7) is a translation with components (x\, x-z). Once 
we choose an orientation on T the complex length is uniquely defined up to sign; 
changing sign corresponds to composing $ with a 180 degree rotation of H 3 taking 
the cylinder to itself. In the limiting case as R — > 00, we obtain a horospherical 
torus; then we define £(7) = for all 7 G ^\{T). 

The Euclidean length L = £(7) of a closed curve 7 on the Euclidean torus T is 
the length of a geodesic in its homotopy class. This is the length of the translation 
M7) = ( x i> x 2), so we have L 2 — x\ + x 2 ., and this can be written terms of the 
complex length as follows: 

(27) L 2 = (coshi? Re£) 2 + (shxhR Im£) 2 . 

Next we extend some geometric notions from simple closed curves on T = Tr to 
arbitrary homology classes in ifi(T;R). 

The complex length of closed curves on T gives a Z-linear function C : H\ (T; Z) = 
7i"i (T) — > C. We define the complex length £ for each element of Hi(T; R) by 
extending this to an R-linear function £ : ffi(T;R) — > C. 

We can regard the generalized Dehn surgery coefficient on T as the homology 
class c = £ _1 (27ri) in Hi(T;M) whenever £ is invertible. (Note that (|26p implies 
that £ is invertible whenever T has tube radius R < 00.) After choosing a basis o, 6 
for Hi(T;Z) = 1? this corresponds to the element (x,y) in R 2 such that xC(a) + 
yC(b) = 2iri, giving the generalized Dehn surgery coefficient as defined in [T5] , 

If dM consists of k tori T\ , . . . Tfc, we have a complex length function Cj : 
Hi(Tj\M.) — * C for each j; the direct sum of these gives a function £ : H\{dM\ R) = 
®jH\(Tj\ R) — > C fc . We then define the generalized Dehn surgery coefficient to be 
the homology class c = £^ 1 (27ri, . . .,2m) G #i(dM;R) = ®jiii (T,; R) whenever 
£ is invertible. 

The Euclidean length L of closed curves on T gives a quadratic form L on 
iii (T; Z) , and this extends naturally to a positive definite quadratic form on iii (T; R} 
If we choose a basis a, b for H\(T\ Z) and let c = pa + qb where p,q G Z, then we 
can write L(pa + g6) 2 = Ap 2 + 2Sp(/ + Cq 2 where A, B, C are constants. We use 
the same formula to define the length L(pa + qb) whenever e I. Then the 
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relationship (|27|) between Euclidean length and complex length continues to hold 
for all homology classes in Hi(T;M,). 

Next we discuss some other geometric quantities that will be important in our 
arguments. 

The area of the torus T = Tr can also be expressed in terms of complex lengths. 
If a, b is a basis for Hi (T; Z), then the area of T is just the area of the parallelogram 
with sides given by the Euclidean translations h(a), h(b). If l a + i9 a and l b + i9 b 
denote the complex lengths of a and b, then 

h(a) = (9 a smhR,l a coshR) and h(b) = (9 b sinhi?, l b coshi?) 

so 



area(Tfl) 



9 a sinh R l a cosh R 
9b sinh R l b cosh R 



sinhi? cosh R (lb9 a — l a 0b), 



provided a, b are oriented so the above determinant is positive. (Note that this is 
the same for any basis a, b which is oriented compatibly with T.) 
We now define the visual area of the torus Tr to be 

(28) A = area(Tfl) 



sinh R cosh R 



This represents the measure of the set of geodesies meeting Tr orthogonally as 
viewed from the core geodesic of the corresponding cylinder of radius R, and will 
play an important role in our analysis. 

Note that the visual area A is the same for parallel tori, i.e. the right hand side 
of (|28jl is independent of R. For a hyperbolic cone manifold with core geodesic of 
length £ and cone angle a we have A = at. In general, A can be expressed in terms 
of the complex lengths on T. Using the notation above we have 

(29) A = l b 9 a - l a 9 b 

for any positively oriented basis a, b for Hi(T; Z). 

Next, we define the normalized length of a homology class c £ Hi(T;M) on 
T = T R to be 

(30) L[c) - L{C) 



y/&re&(T R ) 

This is just the Euclidean length of c after the torus Tr is rescaled to have area 
1. In the case where R = oo, T is a horospherical torus and L is the same for all 
parallel tori, i.e. independent of the choice of horospherical cusp cross section. 

Finally, let M be a hyperbolic 3-manifold with tubular boundary, and let M 
denote its canonical filling. For any sufficiently small r > 0, we can truncate all 
the ends of M to give a hyperbolic manifold with tubular boundary consisting of 
disjoint, embedded tubular tori of tube radius r. Then we define the tube radius R 
of M to be the supremum of all such r. 

Note that if M has boundary components with tube radii R± , . . . , Rk then R is 
larger than R = min(i?i, . . . , Rk)- Further, for any r < R we can truncate M to 
obtain a manifold M with tubular boundary such that all boundary components 
have tube radius r. 
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5.2. Infinitesimal Rigidity keeping Dehn surgery coefficients fixed. 

Now assume we have a hyperbolic 3-manifold M with tubular boundary such that 
each boundary component has tube radius at least Rq = arctanh(l/\/3) ~ 0.65848. 
Then we can use our harmonic deformation theory to prove an infinitesimal rigidity 
theorem for nearby hyperbolic structures: 

Theorem 5.1. Let M be a compact, orientable hyperbolic 3-manifold with tubu- 
lar boundary such that each boundary component has tube radius at least Ro — 
arctanh(l/\/3) ~ 0.65848. Then there are no infinitesimal deformations of this 
hyperbolic structure keeping the Dehn Surgery coefficients constant. 

Proof. Each infinitesimal deformation of the holonomy of a boundary torus T is 
represented, in a neighborhood of T, by a standard form loq with real part 770. We 
can write Wo as a linear combination 

(31) ujq = suj m + (x + iy)uji s,i,!/el 

of the forms u m and uji given in (2) and (3) of [IT]. (The vector fields e2,e3 in [TT] 
are chosen in the directions of the principal curvatures on the tubular boundary.) 
Then the effect of luq on complex length C of any closed peripheral curve is given 
in [TTJ Lemma 2.1] by: 

(32) = ~2sC + 2(x + iy)Re(£) 
at 

and this formula extends by linearity to give the variation in the complex length of 
any element of #i(T;R). 

If the Dehn surgery coefficient c e H\(T; K) is fixed, then the complex length of 
c is C = 27rz. Hence 




= -2s£ = 



and s = 0, so the uj m term vanishes. The contribution to the boundary term 
bnivoi Wo) defined in (|6|) from the boundary torus T was computed explicitly in |11[ 
p382]. Here this simplifies to 

, n 9n sinh R ( 1 \ 

- (x 2 + y 2 ) : — 2 H 5 area(T) < 0, 

v y y coshi?V cosh 2 (i?)y V ; - 

where R is the tube radius of T. 

If M has multiple boundary components, then for any infinitesimal deformation 
keeping the Dehn Surgery coefficients constant, &i?(r/o, ?7o) is a sum of terms of 
this form so &fl(??o,?7o) < 0. But for any non-trivial infinitesimal deformation, 
bnivoi Wo) > by Corollary 14. 31 and we conclude that the infinitesimal deformation 
is trivial. □ 

5.3. Local parametrization by Dehn surgery coefficients. 

Let M be a compact, orientable hyperbolic 3-manifold with tubular boundary, and 
suppose that the tube radius of each boundary component is finite and at least Rq = 
arctanh(l/V3) ~ 0.65848. Let 1Z denote the character variety of representations 
tti(M) — > PS , L 2 (C) up to conjugacy (see [8], [3])- First we describe the local 
structure of the algebraic variety 1Z near the holonomy representation pa : tti (M) — » 
PSL2(C) of M. (Throughout this section we abuse notation, by using the same 
symbol for a representation and its image in the character variety.) 
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Theorem 5.2. Let po : m(M) — ► PSL2(C) be the holonomy representation for a 
compact, orientable hyperbolic 3-manifold M with tubular boundary such that the 
tube radius of each boundary component is finite and at least Rq — arctanh(l/v / 3)- 
Then the character variety 7Z is a smooth manifold near po, of complex dimension 
equal to the number of boundary components ofdM. Further, there is a smooth local 
parametrization oflZ near po by the complex lengths of the Dehn surgery coefficient 
co G Hi (dM ; R) corresponding to po . 

Proof. For each representation p near p , the complex length of peripheral curves 
extends to a well-defined R-linear map 

C p : Hi{dM;R) -► C k , 

where k is the number of tori in dM. Define F : 7Z — > C fc by taking the complex 
lengths of the (initial) Dehn surgery coefficient cq E H±(dM;M.) for M: 

F(p)=£ p (co). 

As in [10;, section 4], the infinitesimal rigidity theorem (Theorem 15. f I) implies that 
the derivative dF po : T po 7Z — > C fe has trivial kernel. Hence the Zariski tangent 
space T po lZ has complex dimension dmicTp 7Z < k. However, by [18\ Theorem 
5.6], we know that 1Z has complex dimension > k at po, so dime T Po 1Z > k. Hence 
dime T po lZ = k and 1Z is a smooth manifold near po of complex dimension k. 
Further, the inverse function theorem then implies that F is a local diffeomorphism. 

□ 

Next, we show that the Dehn surgery coefficients give a smooth local parametriza- 
tion near po. 

Theorem 5.3. Let p Q : ni(M) — > PS , L 2 (C) be the holonomy representation for a 
compact hyperbolic 3-manifold M with tubular boundary such that the tube radius of 
each boundary component is finite and at least Rq = arctanh(l/-\/3). Then there is 
an open neighborhood U of po in 1Z such that for each p € U there is a well defined 
Dehn surgery coefficient c(p) € Hi(dM;M), and the map c : U — ► Hi(dM\R) is a 
diffeomorphism onto its image. 

Proof. For each representation p in a neighborhood V of po in TZ, complex lengths 
define an K-linear map 

£ p : Hi(dM; R) — > 

where k is the number of tori in dM . We can regard this as a function of two 
variables: 

L : V x Hi(dM; R) -> C fc , where £{p, c) = C p (c). 

For each p near po, the corresponding Dehn surgery coefficient c is defined by 
the equation 

(33) C(p,c) = (2m,...,2>Ki). 

For p — po this has a unique solution cq. Differentiating the equation (|33[) at (po, cq) 
gives the linearized equation satisfied by tangent vectors (p, c) to the solution space 
of©: 



dC 
dp 



. dC 

P+jr 

(po,co) oc 



c = 0. 

(po.co) 
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P- 

(po,co) 



Since C(p, c) is a linear function of c this gives: 

dC 
dp 

Now the right hand side is ~dF po (p) where F : 7Z — > C fc is defined by -F(p) = £p(co) 
as in the proof of Theorem 15.21 So this can be written 

£p (c) = -dF P0 (p). 

Since £ Po is invertible and dF Po is invertible by Theorem \5.2\ this has a unique 
solution c for any p and the map p h- > c is invertible. Hence, by the implicit 
function theorem, Q33|) has a unique solution c = c(p) for any p near po and the 
map p > c(p) is local diffeomorphism. Thus the Dehn surgery coordinates give a 
smooth local parametrization of 1Z near po- □ 

This completes the proof of Theorem 1 1.41 

Remark: A hyperbolic structure with infinite tube radius Rj — oo along some 
components Tj of dM can be filled in to give a hyperbolic structure with complete 
cusps corresponding to these components. At such structures, the corresponding 
complex lengths Cj are zero and the corresponding Dehn surgery coefficient is 
defined to be Cj — oo. In this case, a sign for Cj cannot be chosen to vary con- 
tinuously for nearby structures. However, the results of Theorem 15.21 and Theorem 
15.31 extend to this situation provided the complex length is regarded as a function 



C : ® j Hi(T j ;R) -> (C/ ± l) fe , and the Dehn surgery coefficient as an clement of 
® j H 1 (T j ,R) where H^T^.R) = (Hi(Tj,R) Uoo)/±l. 

5.4. An effective version of the hyperbolic Dehn surgery theorem. 

Consider a complete, finite volume, orientable hyperbolic 3-manifold with cusps, 
diffeomorphic to the interior of a compact 3-manifold X with boundary consisting 
of k tori d\X, . . . , dkX . Given a homology class c = (ci , . . . , c^) S H\ (dX; R) = 
Q>jHi(djX;M.), we consider deformations of the hyperbolic structure with Dehn 
surgery coefficients varying "radially" : 

(34) — c, < a < 2tt, 

a 

where a is a smooth increasing function of a parameter t. Then the complex length 
of the homology class Cj for a given value of t is 

(35) C(cj) = a{t)i. 

We want to show that we can deform the hyperbolic structure and increase a 
to 27r, provided the normalized lengths of the surgery coefficients Cj are sufficiently 
large. By Thurston's original Dehn surgery theorem [TH], we can always increase a 
from (corresponding to the complete hyperbolic structure on the interior of X) 
to some small positive value. The local parametrization in Theorem 15 .31 shows that 
we can always increase a slightly, so the set of attainable a is an open subset of 
(0, 27r]. Then we need to control the change in geometry during the deformation, 
and guarantee that no degeneration of hyperbolic structures occurs before a = 2ir 
is reached. 

Here is a brief outline of the argument in this section. Let M t denote the filled 
hyperbolic manifold corresponding to parameter t. By removing disjoint open tubes 
around the ends of Mt we obtain a smooth family of a hyperbolic manifolds Mt with 
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tubular boundary. First we use the positivity condition in Corollary 14.31 provided 
by our harmonic deformation theory to control the variation in complex lengths of 
curves on dM t . This leads to Proposition ^. 51 which gives differential inequalities on 
the total visual area A of the boundary dM t (i.e. the sum of visual areas of all the 
boundary components). Then, in Theorem 15. 6i we apply tube packing arguments 
to obtain a crucial estimate relating A to the tube radius R of M t (as defined at the 
end of section l5~Tj) . This shows that good control on A throughout a deformation 
will guarantee that the tube radius R stays bounded away from zero. Integrating 
the differential inequalities for A shows that such control can be obtained provided 
that the normalized length of the surgery coefficient is sufficiently large (Theorem 
I5.7|) . In Theorem l5.8l we use geometric limit arguments to show that this control on 
the tube radius, together with bounds on the volume (Lemma 15. 9[) and injectivity 
radius of the boundary (Lcmma l5.10j) imply that no degeneration of the hyperbolic 
manifolds M t can occur before a = 2tt is reached. Finally, combining Theorems 15. 71 
and 15.81 gives the main results: Theorems 11.11 and 1 1 . 21 of the introduction. 

Suppose the hyperbolic manifold M t has tubular boundary consisting of tori 
Ti , . . . , Tfe with tube radii R\ , . . . , Rk . We first choose a harmonic representative 
lu for the infinitesimal deformation as in Theorem 14.21 Then near each boundary 
torus Tj we can write the standard part of uj in the form: 

(36) loq = SjU m + (xj + iyj)ui (sj,Xj,yj G M), 
or 

(37) too = Sj{uj m + (Xj -\-iYj)u)i) 
where 

(38) x j =X j s j , a, >,.s,. 

The coefficients Sj , Xj , Yj completely describe the variation in the holonomy of 
the torus Tj for any "radial" deformation of Dchn surgery coefficients as in (f3~4"| . 
In particular, using (133)) and (f3"2"| with C — C(cj), we see that 

1 da 

39 Sj = 

y ' 3 2a dt 

so all Sj are all equal and depend only on the logarithmic derivative of a with 
respect to t. 

Now we choose a parametrization where a is an increasing function of t, with 

Sj = s = — < for all j. Using the crucial positivity property in Corollary 

n-m 2a at 

14.31 we obtain estimates on the size of the coefficients Xj,Yj. In particular, we 
have: 

Proposition 5.4. Let Aj be the visual area of Tj and let A = J2j A? ^ e ^ e 

visual area of dM t . Then if R = min (Ri, ...,R k )>Ro = arctanh(l/ v / 3), 

(40) j2^( x j+o 2 <™ 2 , 

3 

where 

1 , 2 cosh 2 R 

(41) E = t t and w = n n . 

sinh 2 i? (2 cosh 2 i? + 1) sinh 2 i? (2 cosh 2 i? + 1) 
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Proof. By Corollarv l4.31 770 = Rewo satisfies the positivity property 

< b R (r] , 770) = / A ?7o- 

JdM 

Now this integral breaks up into a sum of integrals over the boundary tori Tj, so 
< b R (r] , Vo) = / * D Vo A Vo- 

A priori, some of these boundary integrals could be negative, and this makes the 
argument more complicated in the case of multiple boundary components. 

Note that we have some flexibility in the choice of the tube radii Rj . In particular, 
by adding collars on boundary components we can decrease any Rj, so we may 
assume that Rj = R for each j. By explicit calculations as in p383] we then 
obtain: 

j 

where Aj is the visual area of T.,, s = and 

3 3 2a dt 

-smh 2 R, , o . -2 2cosh 2 i?-l 

(42) a = 2 — (2 cosh 2 i? + l) , b - 



cosh^i? v cosh 2 i?' sinh 2 i?cosh 2 i? 



By the completing the squares we obtain 

< b R ( Vo , Vo ) = Y / a((x j + ±) 2 +YA Ajs 2 + ) A 3 s 2 . 



2a J 1 I J V 4a 

j 

Since a < this gives 

, , \ f Aac~b 2 \ 9 . / 4coshi? \ , , 

43) M%,% < ; ) s A = o 5 s 2 A, 

{ ' W / ; " V 4 « / Vsinh 3 i?(2cosh 2 i?+l)y 

and 

1 a a ^ . b b 2 — Aac 2 

where we write A — > „■ Ai . Computing — = £ and ^ — = w 

^ J J 2a 4a 2 

gives the result. □ 



Combining this result with equation (|32() gives us control on the holonomy of 
peripheral curves: the variation in the complex length C of any homology class on 
a boundary torus Tj is given by 

(45) ^-(C) = - 2sC + 2 ( x j + Wo) Re (£) = -^r-£ + 2(afj + iy*) Re(£). 

a ai 

We now use this to estimate the variation in the visual area Aj of Tj. Choose 
an oriented basis a, 6 for Hi(Tj)Z) with complex lengths l a + i9 a ,lb + i&b- Then 
applying the formula (|45|) to a, 6 gives 

din dd a 1 da ,, 
46 -JL+i-JL = - — {l a + i6 a ) + 2{ Xj + i V j)l a 

dt dt a dt 
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and 

dlh dOh Ida,, 
(47) ^ +i -A = __ {lb + l b)+2{xj+Wj)l , 

By differentiating equation it follows that 

(48 » |(^) = |(w--w = ^Qf + ! 



< 49 > 

where the X,- satisfy the inequality (|40D . 

In the following argument, the iotaZ visual area of the boundary A = ■ »4j will 
play a crucial role. To control the behavior of A, we will consider the variation of 

(50) «i = = ^ = I>i- 

3 

Note that these quantities only depend on the canonical filling M t ; any truncation 
M t of M t with tubular boundary gives the same values for Aj, A, Vj and v. 

First we examine the limiting behavior of Vj for our family of hyperbolic man- 
ifolds with tubular boundary M t , where we fix a homology class Cj € H\ (Tj ; R) 
and vary the Dehn surgery coefficients — Cj along a ray going out to oo as a — > 0. 
The homology class Cj has complex length C(cj) — ai where a is the deforma- 
tion parameter. Hence, by equation (P7)l . the Euclidean length of Cj on Tj is 
Lj = L(cj) = asmhRj where Rj is the tube radius of Tj, so 



area(Tj ) Aj sinh Rj cosh R 



Lj a 2 sinh (Rj ) 



Vj coth Rj 



As a — > 0, the hyperbolic structures converge to the complete hyperbolic structure 



on the interior of X and Rj — * oo. Hence 



area(Tj ) 



tanh R.j — > ^— . 



and 

3 3 L 3 

where Lj — L(cj) is the normalized Euclidean length of Cj, as defined in on a 
horospherical cross section for the complete hyperbolic structure on the interior of 
X. 

Remark: The quantity L = L(c) defined by 

(5i) i = yi 

K ' L 2 L 2 

3 3 

seems to be a useful analogue of the normalized length in the one-cusped case, and 
we will also call it the normalized length of the homology class c S Hi(dX;M) in 
the multi-cusped case. Its reciprocal 1/L gives a good measure of the distance from 
the complete hyperbolic structure on the interior of X to the hyperbolic structure 
M(c) with Dehn surgery coefficient c. 
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Differentiating vj using equation (|48|) gives 

dvj _ d ( Aj\ _ 1 ( dAj 2Ajda\ _ 2AjXj 

{5) ~dT ~ Tt\tf)~tf\dt ~a~~dt)~~tf~~ jXj 

and 

(53) i^ = 2 y v -i Xj = -l^y v -i Xj . 

y ' v dt ^ v 3 a dt ^ v 3 

j j 

Combining this with Proposition 15.41 gives our basic differential inequalities for v. 

Proposition 5.5. Let v = A/ a 2 = v± + . . .+Vk where Vj = Aj/a 2 . Then v satisfies 
the differential inequalities 

1 Ida ldv 1 /2cosh 2 7?-l\ Ida 

(54) > 7.-JI > 



sinh 2 R a dt v dt sinh 2 7? ^ 2 cosh 2 _R + IJ a dt 1 

provided the tube radius R of the canonical filling M t is larger than Rq . Further, v 
satisfies the initial condition 

(55) lim v — -z—, 

3 3 

where Lj = L(cj) is the normalized Euclidean length of cj, on a horo spherical cross 
section for the complete hyperbolic structure on the interior of X. 

Proof. Since R > Rq, we can truncate M t to give a hyperbolic manifold with 
tubular boundary M t such that all components of dM t have tube radius at least 
R, where R> R> R . 

Using the Cauchy-Schwartz inequality and (|40|) then gives 



v z — ' v z — ' v 

3 3 3 



since J^j v j = v - Hence, 

(56) -w-£<^2 ~ x 3 < W -C 



v 

3 



Writing out w and £ in terms of R gives 



1 v-^ Vj 1 /2cosh 2 i?-l 



sinh 2 i? 4^ v 3 sinh 2 i? \2cosh 2 J R+ 1 

3 

Multiplying through by the negative number — and recalling, from (|53[) , that 

a dt 



1 dv 1 da 
v dt a dt '— J v 



E-v. 

gives the inequality 



1 1 da > ldv > 1 /2cosh 2 i?- 1\ 1 <!n 



sinh R a dt v dt sinh R V 2 cosh R +1/ a dt 
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Since these inequalities hold for all R such that Rq < R < R, they also hold when 
R is replaced by R. This gives (|54l) . The initial conditions for v was already derived 
above. □ 

Next, let M be a hyperbolic 3-manifold with tubular boundary, and let M denote 
its canonical filling. Then the tube packing arguments in the proof of [111 Theorem 
4.4] give us the following crucial estimate relating the total visual area of dM to 
the tube radius of M . 

Theorem 5.6. Let M be a compact, orientable hyperbolic 3-manifold with tubular 
boundary, and let A be the total visual area of dM. Then 

(57) A>h(R), 

where R is the tube radius of the filled manifold M , and h is the function given by 

(58) Mr) = 3.3957^^-. 
v ' w cosh(2r) 

Proof. We briefly recall the argument from [TT]. If we expand tubes around the 
ends of M at the same rate, then these first bump when the tube radius is R. 
Suppose that the tube bounded by a torus Ti bumps into the tube bounded by 
a torus Tj when the tube radius reaches R; possibly with i = j. Then the tube 
packing arguments from [11] show that Ti U Tj contains two open ellipses meeting 
only at the bumping point, each with semi-major axes 

0.980258 sinh R cosh R , , sinh R sinh R 



a = 



and b = 

cosh(2i?) sinh(2i?) 



and hence of area 



Af, = nab 



0.980258 7T sinh 2 R 



2cosh(2i?) 

(These two ellipses are in the same torus Tj if j = i; otherwise there is one ellipse in 
Ti and one ellipse in Tj.) Further, the packing density for these ellipses is at most 
so it follows that 



2V3' 



. 4^ , sinh 2 i? 
areaffi U Tj) > -^A e > 3.3957 — 

V " ~ TT ~ COSh(2i?) 



and the total visual area of dM satisfies 



A > > 3.3957 , = h(R). 

sinh R cosh R cosh R cosh(2i?) 



□ 



Now we can apply the same arguments as in [ITJ section 5], but with the tube 
radius condition R > 0.531 replaced by R > Rq = arctanh(l/ v3) ~ 0.65848. Since 
h(r) is a decreasing function for r > Ro it follows from (|57|l that if initially the tube 
radius satisfies R > Rq and we know that A < h(Ro) throughout a deformation, 
then R > Rq throughout the deformation. 

Next we use the control on A given by Proposition l5.5l and the inequality (|57p to 
show that the tube radius R of M t stays bounded below throughout any deformation 
as in (f3~4"|) with < a(t) < ao < 2tt, provided that the normalized lengths of the 
Dehn surgery coefficients Cj G Hi(Tj\M.) are sufficiently large. 
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First note that the inequalities (|54| are exactly equivalent to the differential 
inequalities for u — 1/v obtained in |f f j in the one cusped case: 

f 1 da Idu 1 ( 2 cosh 2 ^ - 1 \ Ida 

< < 



sinh 2 R a dt u dt sinh 2 fl I 2 cosh 2 E + 1 J a dt ' 
Now we choose a parametrization with t — a 2 . Then this becomes 
1 u du 1 /2cosh 2 i?— l\ u 



( ' )J1 sinh 2 R 2a 2 ~ dt ~ sinh 2 R \ 2 cosh 2 R + 1 J 2a 2 ' 

We analyze this as in |f f j by introducing a new variable z = tanh(p) where 
h(p) = A and p > Rq. Note that p is defined whenever A < h(Ro) and if R > Ro 
then (|57| implies that Ro < p < R- This allows us to replace R by p in the 
inequality ([59]) . 

Now we define functions 

(60) = A- = W) = 3.3957z(l-z 2 ) ' 

r«-n a H(z)i-#_ l + z 2 

(61) —=6J9il^> 

and 

(62) ew = F W (l-^)(l +2 »)_ (l + z 2 ) 2 



2 z 2 (3-z 2 ) 6.79f4 z 3 (3-z 2 )' 

Then the differential inequality (|59[) for m, with i? replaced by p, becomes 

-ow<f saw 

dz 

and putting u = tH(z) gives the differential inequalities for — obtained in [ITJ 
equation (51)]: 

g(g) rfz < 1 < H'(z) dz 
H{z) - G(z) dt ~ t ~ H(z) + G{z) dt ' 

To solve these differential inequalities we write: 

H'(z) ~, , 1 H'(z) , . 1 

- F(z) + , - , ' r = F(z 1 



H(z)-G{z) l-z 1 H(z) + G(z) 1-z 

Then 

_ (l + 4z + 6z 2 + z 4 ) 
(z + l)(l + z 2 ) 2 
is integrable on the interval < z < 1 and 

z 6 + 7z 4 + 12z 3 - 9z 2 - 4z + 1 



F{z) 



(z + 1) (z 2 + 1) (z 2 -22-1) (z 2 + 2z - 1) 



is integrable on the interval \[2 — 1 + e < z < 1 for each e > 0. 

Now we integrate (|53")) with respect to i from t\ to t where < t\ < t < 1, 
and z(r) > tanh(i?o) = 1/V3. Carefully taking a limit as ii — > 0, zi = z(£i) — > 1 
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as in gives estimates on the time t — a 2 taken to reach a given value of 
A = A(z) = jjjjj where z = z(r). Define L by 

where Lj = L(cj) is the normalized length of the homology class Cj £ Hi(Tj-,M.) for 
the complete hyperbolic structure on the interior of X, as in (|5ip . Then we obtain 



(64) 

where 

(65) 

and 

(66) 



/(*) > j- > /(*), 



f{z) = 3.3957(1 - z) exp(- J F{w) dw) 



f(z) = 3.3957(1 - z) exp(- J F{w)dw). 



These bounds are illustrated in the graph in Figure 1. 




0.2 0.4 0.6 0.8 1 1.2 1.4 



Figure 1 . Graph of A versus x = 



We conclude that we can increase the parameter a from to 2-7T, maintaining 
z = tanh p > zq = tanh(i?o), hence keeping the tube radius R > p > Ro = 
arctanh(l/v / 3) and A < h(Ro), provided 



L > 



^ exp f f ° F(w)dw I w 57.5041 

3.3957(1 -z ) y \J 1 V ' 1 



or 



L > V57.5041 w 7.58315. 



Thus, we have shown that as long as L satisfies this inequality then there is a lower 
bound to the tube radius: 

Theorem 5.7. Let M t be a smooth family of hyperbolic structures with tubular 
boundary on X for < t < 1, with Dehn surgery coefficients ^fyC where c = 
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(ci, . . . , Cfe) 6 Hi(dX;M) and < a(t) < ao < 2ir. If the normalized lengths of the 
surgery coefficients Li = L(ci) satisfy 

V J- < 4? where C = 7.5832, 

i i 

then tube radius of the filled manifold M t satisfies R > Rq for all t. 

The next main result is the following analogue of [11] Theorem 5.4]: 

Theorem 5.8. Let M t be a smooth family of hyperbolic structures with tubular 
boundary on X for < t < 1, with Dehn surgery coefficients ^fyC where c = 
(ci, . . . ,Cfe) S Hi(dX;M) and a(t) is an increasing function oft with < a(t) < 
a < 2t:. Suppose the tube radius of the canonical filling M t satisfies R > Rq = 
arctanh(l/v / 3) ~ 0.6585 for t = and the total visual area of dM t satisfies A < 
h{) = h(Ro) for all t. Then the manifolds M t converge geometrically as t — * 1 
to a hyperbolic manifold M% with tubular boundary. Further, the Dehn surgery 
coefficient for Mt converges to the Dehn surgery coefficient for Mi . 

Proof. We begin with some estimates on the geometry of the manifolds M t . First 
we study the change in volume of the filled in manifolds M t with Dehn surgery type 
singularities. In fact, we have Schlafli type formula for the variation in volume: 

Lemma 5.9. Let M t be a smooth family of hyperbolic 3-manifolds with tubular 
boundary with Dehn surgery coefficients varying radially as in {3$ . Then the vari- 
ation in volume V of the filled manifolds M t is given by 

A 

(67) dV = -—da, 

2a 

where A = ■ Aj is the total visual area of dM t . In particular, the volume de- 
creases as a increases. 

Remark: In the cone manifold case, (l67|) is just the usual Schlafli formula: 

dV = —\'Yli-ijda, where £j is the length and a is the cone angle of the core 

geodesic produced when Tj is filled. 

Proof. From [T^l equation (46)] (or chapter 5], pj section 4.5]) the variation in 
volume is a sum of contributions dVj from the boundary tori Tj, and we have 

(68) dVj = -^{l b d8 a -l a d9 b ), 

if a, b is any oriented basis for Hi(Tj;Z). Now, from equations (14^)) and (|4T|) . we 
have 

d6 a _ d£i 
b dt a dt 



Hence 

dV 

□ 



, . 1 da „ „ , . , . 1 da . „ , , 
h(--j-0 a + 2yl a ) - l a (--j-0 b + 2yl b ) 
a at a at 

1 da (l ft I \ - A i da 
a dt a dt 



E^-da = — — da. 
2a 2a 
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Tube packing arguments as in give the following estimate on injectivity 
radius of the boundary of a hyperbolic 3-manifold with tubular boundary. 

Lemma 5.10. Let M be a compact, orientable hyperbolic 3-manifold with tubular 
boundary such that each boundary component has tube radius at least R > 0. Then 
there exists a constant c(R) > such that M can be truncated along tubular tori to 
give a hyperbolic 3-manifold with tubular boundary such that the injectivity radius 
of the Euclidean metric on each boundary component is at least c(R) . 

Proof. We first expand the boundary tori of dM , moving each torus inwards at the 
same rate until it bumps into itself or another boundary component. Let T\, . . . , T k 
be the (immersed) tubular tori obtained when this bumping occurs; these tori meet 
tangentially in a finite collection of points. Let Ri be the tube radius of X^; then 
Ri> R for all i. 

We let N — M denote the canonical filling of M, N the universal covering of N, 
and N the metric completion of N. For each subset Y C N C N, let Y denote the 
closure of Y in N. 

We can regard Ti, . . . , T k as a subsets of the canonical filling N. Then each Tj 
bounds an open tube V, in N and the tubes Vy, . . . , Vk are disjoint. 

Let Ui be one (fixed) lift of Vi to the universal cover N of N and let f/j be its 
closure in N. Now consider the lifts U of tubes Vj to N such that dU meets dUi 
tangentially at a point (possibly with j = i), and let U be the closure of U in N. For 
each such U, we construct a point q S U as follows: Let p denote the intersection 
point of dUi and dU . Then q is the point inside U at distance R from p along the 
geodesic through p orthogonal to dU. 

Claim. Let Q be the collection of all points q constructed as above. Then the 
distance in N between any two distinct points of Q is at least 2R. 

Proof. Let V, U" C N be lifts of the tubes Vi,...,V k whose closures U', U" C N 
contain two distinct points q',q" G Q, and let B(q',R) and B(q",R) denote the 
open balls in N of radius R around q' and q" respectively. Since R < Rj, it follows 
from the triangle inequality that B(q',R) C U' and B{q",R) C U". But U 1 and 
U" have disjoint interiors, hence B(q',R) and B(q",R) are disjoint. This proves 
the claim. □ 

Now the distance from the boundary of Vi to any singular point in the completion 
of N is at least R. So we can expand Ui to an open tube Wi of radius di — Ri + R 
which embeds isometrically in N. 

The geometry of Wi can be described as follows. Let g be a geodesic in H 3 . The 
universal cover of H 3 — g can be completed by adding a geodesic, g, which projects 
to g in H 3 . (This can be thought of as the infinite cyclic branched cover of H 3 
branched over the geodesic g.) Let H 3 denote this completion. Then there is an 
isometry (f) : Wi — > W to the open tube W of radius di about g in H 3 . 

We use this to identify Wi with W C H 3 and to identify the set of points Q with 
a subset Q C dW C H 3 . Since the closure W of W is a convex subset of H 3 , it 
follows from the claim above that d^ 3 (q' , q") = dyy(q', q") > 2R for all q' ^ q" E Q. 

Next we use the arguments of [TTj to estimate the distance measured on dUi 
between the tangency points p described above. For q S Q, let B q denote the 
closed ball in H 3 of radius R around the point q and let P q denote the orthogonal 
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projection of B q onto the surface at radius Ri from the core geodesic of H 3 . Since 
the balls B q have disjoint interiors, are of equal radius and are allcentered at the 
same distance di — Ri + R from the core geodesic in H 3 , it follows that their 
projections P q also have disjoint interiors. 

As in [11] we introduce cylindrical coordinates on H 3 around the geodesic g and 
lift these to cylindrical coordinates (r, 9, £) on H 3 ; here the angle is a well defined 
real number. Then from [ITJ Lemma 4.3], the projection of B q to the (6, Q plane 
can be parametrized by 

sinh 2 C cosh 2 (i?i + R) + sin 2 6 sinh 2 (i? 4 + R) < sinh 2 i?, 

where R < Ri. 

Now, as in [IT], we have 

, , sinh R sinh R 1 

| sinhCl < — < — , , om < 



cosh(Ri + R) ~ cosh(2i?) ~ 2\/2' 
hence, using convexity of the sinh function, 

IsinhCl < S\C\ where S = 2V2 = < . 

1 Sl " ISI arcsinh-^ 0.980258 

2V2 

Further, 

I sin 6»| < \B\, 

so P q contains the region: 

(SQ 2 cosh 2 ^ +R) + 6 2 sinh 2 (i? l + R) < sinh 2 i? 



S cosh(i?, + R) \ 2 , _ , x2 f smh(Ri+R)\ 

1 (Ccoshi?,) 2 + v ; (9 sinh Ri) 2 < 1. 



sinh R cosh Ri J \ sinh R sinh Ri 

Since ( cosh Ri and sinh i?^ are Euclidean coordinates on the surface in H 3 at 
radius r = Ri, this equation describes an ellipse with semi-major axes 

sinh R cosh Ri sinh R sinh Ri 

and 



S cosh(i?i + R) smk{Ri + R) 

whose interior is disjoint from the ellipses corresponding to other points q G Q. 
Thus these project to ellipses in Ti with disjoint interiors. 
Now 

- = S(cothi? + tanhi? 4 ) < S"(cothi?+ 1) 
a 

and 

- = coth R + coth R t < 2cothi? < cothi?+ 1 < 5(cothi?+ 1), 


where S < 980258 ' Hence the Euclidean injectivity radius of is larger than 

0.980258 



c(R) 



coth R+l 



By adding sufficiently small collars on the tori Ti we obtain a hyperbolic manifold 
with tubular boundary as desired. □ 
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The proof of theorem 15.81 now follows from the arguments in the proofs of [TT| 
Theorem 1.2] and [TTJ Theorem 3.12] together with the volume and injectivity 
radius estimates given in Lemmas 15.91 and 15. 101 above. We outline the argument. 

Let M be a hyperbolic 3-manifold with tubular boundary whose canonical filling 
M has tube radius at least Rq- Then by using Lemma \5. 101 we can add standard 
collars on the boundary components to obtain a hyperbolic manifold N C M with 
tubular boundary such that the boundary components of dN are at least distance 
Ro/2 apart and have tube radius at least Ro/2, and whose injectivity radius at 
all boundary points (as defined in [TTJ p388]) is at least d, where d is a positive 
constant depending only on Rq. 

Thus we can truncate the manifolds M t to obtain a smooth family of hyperbolic 
manifolds N t with tubular boundary, with uniform lower bounds on the tube radius 
and injectivity radius (as defined in 11, p388]) on their boundaries. Further, the 
volume of Nt is at most the volume of the canonical filling Nt = Mt, an d Lemma l5~9l 
shows that this is decreasing throughout the deformation since a(t) is increasing. 
Hence the volumes of the Nt are bounded above. 

It follows as in the proof of [TTJ Theorem 3.12] that Nt converge in the bilipschitz 
topology to a hyperbolic manifold N± with tubular boundary, and we can choose 
holonomy representations pt for Nt converging to the holonomy representation p\ 
for N\. Further, the convergence of the geometry on the boundary of Nt implies, 
using (|26p . that the corresponding complex length functions Ct for dN t converge 
to the complex length function Ct for dN%. 

In particular, since lim^i £t(cj) = limt— >i a(t)i ^ the limiting complex length 
function is not identically zero. Hence the limiting tube radius for each boundary 
component of Ni is finite. It follows that L\ is invertible and the Dehn surgery 
coefficient for each boundary component of N\ is the limit of the corresponding 
Dehn surgery coefficients for N t . This completes the proof of Theorem 15.81 □ 

Combining the control on tube radii given by Theorem 15.71 with Theorem 15.81 
shows that we can deform the complete hyperbolic structure on the interior of X to 
a hyperbolic structure with Dehn surgery coefficient c E Hi(dX;R.) provided the 
normalized length L of c is at least 7.5832. This proves our main result. 

Theorem 5.11. Consider a complete, finite volume hyperbolic structure on the 
interior of a compact, orientable 3-manifold X with k > 1 torus boundary compo- 
nents. Let Ti, . . . , Tfc be horospherical tori which are embedded as a cross-sections 
to the cusps of the complete structure. Then there exists a universal constant 
C = 7.5832 such that there is a "radial" deformation from the complete hyperbolic 
structure on the interior of X to a hyperbolic structure with Dehn surgery coefficient 
c = (cx,C2,-- - ,Cfc) € Hi(dX;WL) through hyperbolic structures with Dehn surgery 
coefficients ^-c provided the normalized lengths Lj = L(cj) on Tj satisfy 



In particular, this implies Theorems 11.11 and 1 1 . 21 as stated in the introduction. 

Remark: The same result holds if some surgery coefficients Cj are infinite, i.e. 
some cusps remain complete. The proof is essentially the same as given above. 
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5.5. Volume estimates. 

Finally, we use our estimates on the total visual area A to control the change 
in geometry during hyperbolic Dehn filling. We consider a deformation as in the 
previous section, and use the notation and results from that section. 

The inequalities in (|64[) give us upper and lower bounds on A as a function of 
x = provided x < ( 7 lg 32 ) 2 - If L > 7.5832 then these estimates apply throughout 
the deformation, as a increases from to 2-7T. Then by integrating the estimates, 
we obtain upper and lower bounds on the change in the volume V of the filled in 
manifolds Mt with Dehn surgery type singularities. 

Using the parametrization t — a 2 , Schlafli's formula (l67j) gives 

d . , A da A A 
~dV '~~2cVdl"^\^^^\~t 

2 + 

Fixing L and putting x = j-^ = we find that the decrease in volume is 

Jo 4f J a Ax 

where x — ^rj- and A{x) lies in the region defined by the inequalities in (|64[) and 
illustrated in Figure [JJ 

The lower bound for A(x) is given by A — A(z) where x = f(z), dx = f'(z)dz. 
Hence 

1 ? J, <AV 



4 Ji f(z) 

where f(z) = x = ^ 2 . Rewriting this using the definitions of / and A = 1/H 
gives 

1 f 1 H'(z) 



■ dz < AV. 



4 J- z H(z)(H(z)-G(z)) 
The upper bound for A(x) is given by A — A(z) where x = f(z), dx = f'(z)dz. 
Hence 



Av<u z ^m d z 



where f(z) — x — . Hence 



4 Ji f(z 

2 

L 2 



A J, H(z)(H(z) + G(z)) 
In particular, taking £ = l/v3 gives an upper bound AV < 0.197816 when 

^ff- < /(1/V3), i.e. L > 7.5832. 

This gives us the following estimates on the changes in geometry during gener- 
alized hyperbolic Dehn filling. 

Theorem 5.12. Let X be a compact, orientable 3-manifold as in Theorem \5.11l 
and let Vco denote the volume of the complete hyperbolic structure on the interior 
of X . Let c — (ci,...,Cfc) € Hi(dX;M) be a surgery coefficient with normalized 
lengths Li — L{ci) satisfying 

J- = V i- < i where C = 7.5832, 
L 2 L 2 C 2 

3 J 
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and let M(c) be the filled hyperbolic manifold with Dehn surgery coefficient c. Then: 
(1) The decrease in volume in hyperbolic Dehn filling AV — Voo — vol(M(c)) 
satisfies 

1 dz<AV<- — , , dz, 



4Jz H{z)(H{z)-G(z)) ~ - 4/g H{z){H(z) + G(z)) 
where z and z are defined by 

f{z) = f(z) = 

(2) The total visual area A in M(c) satisfies 

A*) = 777^7 <A< A(z) ' 



where z, z are as above. 

Remark: The estimate on the change in volume during Dehn filling in Theorem 
15.121 is a significant improvement on the estimate obtained in 11. The previous 
analysis bounded the change in volume until A = a£ reached its maximum allowed 
value h m axi but in this process the parameter a could increase beyond 27r. Here 
we estimate the change in volume until a reaches 2-7r; this gives the more refined 
estimate. 

The graphs in Figures 2 and 3 illustrate the results in Theorem l5.12l The dotted 
lines in these figures correspond to the asymptotic formulas of Neumann-Zagier [15j : 
as L — > 0Oj the decrease in volume is AV ~ and the visual area is A ~ v . 




0.1 0.2 0.3 0.4 0.5 0.6 



Figure 2. Graph of AV versus x = 

In particular, taking L = 7.5832 in the above gives the following numerical 
estimates. 

Corollary 5.13. Let X be a compact, orientable 3-manifold as in Theorem \5.11\ 
and let c S H±(dX;M.) be a surgery coefficient with L{c) > 7.5832. Then 

(1) the decrease in volume during hyperbolic Dehn filling is at most 0.198, 

(2) the total visual area of the boundary of the filled hyperbolic manifold M(c) 
is at most h(R ) w 0.980254. 

In particular, if M{c) is a hyperbolic manifold with a smooth core geodesic of length 
£ then A = 2n£. Thus the core geodesic length is at most ~ 0.156012. 
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